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ABSTRACT 


- ^This  dissertation  describes  the  development  of  new 

techniques  for  calibrating  a  microwave  automatic  network 
analyzer  (hereafter  ANA)  for  complex  reflection  measure¬ 
ments.  These  techniques  promise  significant  advantages 
for  miniature  connector  systems.  They  include  the  use  of 
an  experimentally-derived  characterization  of  an  APC-7 
open-circuit  termination  as  a  high-reflection  calibration 
standard,  adapter  correction  based  on  low  dissipative 
losses,  and  a  length  of  transmission  line  instead  of  a 
sliding  load  or  fixed  termination  as  a  low-reflection 
standard. 


Gains  afforded  by  jnese  new  techniques  are  reduced 
operator  intervention  during  calibration,  less  wear  of 
the  measurement  port  connectors,  lower  equipment  and 
maintenance  costs,  and  calibrations  referenced  to  standards 
truly  relevant  to  the  interconnecting  transmission  medium. 
Computer  memory  requirements  are  increased,  but  the  cost 
of  the  increase  is  insignificant  compared  to  the  cost  of 
microwave  instrumentation.  Slightly  longer  calibration 
times  are  offset  by  reduced  demands  upon  operator  skill. 


Methods  of  extending  these  calibration  techniques 
to  full  two-port  measurements  will  be  discussed. 
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ABSTRACT 

This  dissertation  describes  the  development  of  new 
techniques  for  calibrating  a  microwave  automatic  network 
analyzer  (hereafter  ANA)  for  complex  reflection  measure¬ 
ments.  These  techniques  promise  significant  advantages 
for  miniature  connector  systems.  They  include  the  use  of 
an  experimentally-derived  characterization  of  an  APC-7 
open-circuit  termination  as  a  high-reflection  calibration 
standard,  adapter  correction  based  on  low  dissipative 
losses,  and  a  length  of  transmission  line  instead  of  a 
sliding  load  or  fixed  termination  as  a  low- re flection 
standard. 

Gains  afforded  by  these  new  techniques  are  reduced 
operator  intervention  during  calibration,  less  wear  of 
the  measurement  port  connectors,  lower  equipment  and 
maintenance  costs,  and  calibrations  referenced  to  standards 
truly  relevant  to  the  interconnecting  transmission  medium. 
Computer  memory  requirements  are  increased,  but  the  cost 
of  the  increase  is  insignificant  compared  to  the  cost  of 
microwave  instrumentation.  Slightly  longer  calibration 
times  are  offset  by  reduced  demands  upon  operator  skill. 

Methods  of  extending  these  calibration  techniques 
to  full  two-port  measurements  will  be  discussed. 
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INTRODUCTION 

The  HP  8542B  ANA  has  been  the  laboratory  standard 
of  one-and  two-port  microwave  measurement  systems  for 
frequencies  up  to  18  GHz  since  its  introduction  to  the 
electronics  marketplace  in  1969.  Its  function  is  to  make 
highly-repeatable  frequency-domain  measurements  of  precise 
standard  terminations ,  to  compute  a  set  of  error  cor¬ 
rection  coefficients  from  these  data,  to  apply  these  cor¬ 
rections  to  measurement  data  of  unknown  devices,  and  to 
output  the  results  in  one  or  more  of  the  many  tabular  or 
graphic  formats  available.  In  spite  of  the  ANA's  twelve- 
year  availability,  its  calibration  procedures  have  re¬ 
mained  essentially  unchanged. 

The  purpose  of  this  dissertation  is  to  present 
innovative  calibration  techniques  developed  during  the 
research  for  this  project.  These  techniques  offer  sig¬ 
nificant  benefits  compared  to  the  previous  calibration 
procedures . 

For  one,  an  experimentally-derived  characterization 
of  an  APC-7  open  circuit's  frequency  dependence,  which 
very  closely  corroborates  a  theoretical  characterization 
(Ref.  1) ,  permits  the  elimination  of  offset  short  circuits 
as  high-reflection  calibration  standards.  Open  circuits 
offer  at  least  two  major  advantages  over  offset  shorts. 
First,  a  single  open  circuit  can  be  used  in  a  wide-band 
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calibration  (i.e.,  2-18  GHz),  obviating  the  need  for 
connection  of  several  offset  shorts  (one  for  each  band, 
2-4,  4-8,  8-12.4,  12.4-18  GHz).  As  a  result,  wear  of  the 
connectors  at  the  measurement  port  can  be  reduced  sig¬ 
nificantly.  Connector  wear  is  the  fundamental  limit  on 
the  usefulness  of  the  "stored  calibration"  principle. 

Second,  the  length  of  each  offset  short  represents 
the  X/4  plane  only  at  the  center  frequency  of  each  band. 
Thus,  the  use  of  offset  shorts  can  lead  to  band-edge  loss 
of  accuracy. 

Further,  the  adoption  and  implementation  of  an 
adapter  correction  procedure  theorized  by  A.  Uhlir,  Jr. 
(Ref.  2) ,  provided  a  means  by  which  an  APC-7  calibration 
can  be  extended  for  use  in  any  other  transmission  format 
without  the  need  for  a  complete  calibration  kit  in  that 
format.  Calibration  kits  are  very  expensive  and,  in  the 
case  of  SMA  (Figs.  90  and  91  show  sectioned  views  of 
typical  SMA  connectors  and  line) ,  the  components  are 
extremely  delicate,  sensitive  to  wear,  and  easily  damaged 
even  when  handled  with  care.  Their  use  also  shortens  the 
useful  life  of  the  adapter  as  a  tool  for  precise  measure¬ 
ments. 


Furthermore,  the  employment  of  a  fixed  load  and  a 
fixed  length  of  transmission  line  in  the  exact  measurement 
medium,  as  a  replacement  for  the  air  sliding  load  as  a 
zero-reflection  standard,  eliminates  the  discontinuity 
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errors  or  ambiguity  inherent  in  attempting  to  define  zero- 
reflection  in  a  dielectric  filled  medium  (i.e.,  SMA)  with 
an  air-dielectric  standard.  In  addition,  the  incorpo¬ 
ration  of  a  double  running  average  technique  to  cancel  the 
effect  of  reflection  residuals  of  the  load  and  line  e- 
liminates  the  errors  (or,  more  rarely,  divergences)  in  the 
circle-fitting  routines  incorporated  in  the  standard 
software  for  the  system. 

The  combined  effect  of  these  developments  is  to 
lessen  operator  interaction,  to  decrease  wear  of  the 
measurement  port  connectors,  to  reduce  the  number  of 
required  calibration  standards,  and  to  diminish  certain 
systematic  errors.  Thus,  the  changes  presented  here  im¬ 
prove  the  functioning  of  automatic  microwave  network  ana¬ 
lyzers  and  extend  their  applicability  to  a  wider  variety 
of  transmission  media. 
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GENERAL  BACKGROUND 

The  Hewlett-Packard  HP  8542B  automatic  microwave 
network  analyzer  (ANA)  is  the  basis  for  this  research 
although  the  techniques  developed  here  can  readily  be 
adapted  for  use  with  other  microwave  measurement  systems. 
This  ANA  is  a  versatile  transfer-function  measurement 
system  which  can  make  rapid ,  highly  accurate,  repeatable 
measurements  of  scattering  parameters  from  110  MHz  to  18 
GHz  in  all  coaxial  and  waveguide  formats  for  which 
standard  impedances  are  available  for  calibration  purposes. 
This  research  broadens  the  domain  of  suitable  standards. 

All  the  research  for  this  document  was  conducted 
using  the  HP  8542A  ANA  at  M.I .T.-Linclon  Laboratory  which 
has  been  updated  and  modified  to  be  essentially  equivalent 
to  an  HP  8542B.  The  most  important  feature  of  this 
updating  is  crystal-controlled  frequency  synthesis. 

In  addition,  this  system  contains  several  options 
which  make  it  as  powerful  and  convenient  as  any  more 
recently  developed  systems.  Three  options  which  proved 
invaluable  throughout  this  project  are  the  8500A  System 
Console  (Maxi  System) ,  the  TODS-II  Test  Oriented  Disc 
System,  and  the  Versatec  Matrix  200  plotter. 
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The  8500A  System  Console  consists  of  a  keyboard, 
control  panel,  interface  unit,  display  generator,  large- 
screen  CRT  display,  track  ball,  and  alpha-numeric  printer. 
The  TODS-II  Disc  System  provides  storage  of  and  ready 
access  to  source  and  data  files.  A  core  resident  monitor 
program  performs  directory  search,  program  loading,  and 
control  of  the  operator  terminal.  A  hard  copy  of  both 
alphanumeric  and  graphic  data  can  be  obtained  from  the 
Versatec  Matrix  200  Plotter.  A  block  diagram  of  the 
total  system  is  displayed  in  Figure  1. 

All  the  controlling  computer  programs  were  written 
in  Hewlett-Packard  combined  ATS  and  ANA  BASIC  which 
offered  full  measurement,  computational  and  output  flexi¬ 
bility,  required  no  compilation,  and  enabled  instantaneous 
line-by-line  editing  from  the  operator  console.  The 
programs  assume  that  the  TODS-II  Disc  System  is  available 
and  require  at  least  8K  of  usable  core  memory.  They 
perform  all  the  necessary  manipulations  for  storing  data, 
controlling  RF  measurements,  computing  correction  coef¬ 
ficients,  applying  the  corrections  to  Device  Under  Test 
(DUT)  data,  and  outputing  the  results  in  an  appropriate 
format  for  meaningful  interpretation. 


'«#».*  ii  IfciMI 


The  use  of  BASIC  will  permit  the  programs  to  be 
adapted  to  newer  ANA  models  which  are  exclusively 
programmed  in  variants  of  BASIC.  (The  HP  8542B  would 
operate  faster  and  have  more  available  memory  if  the 
programs  were  written  in  FORTRAN.) 
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FREQUENCY  CHARACTERIZATION  OF  APC-7  OPEN-CIRCUIT 
REFLECTION 

Open-circuit  coaxial  standards  have  been  generally 
used  only  at  relatively  low  frequencies  (below  2  GHz) 
where  frequency  dependence  of  phase  can  be  approximated 
as  linear.  It  appeared  that  an  open  circuit  would  be  an 
excellent  reflection  standard  at  all  the  usable  frequenc¬ 
ies  of  the  ANA  (up  to  18  GHz)  if  a  highly  accurate, 
reproducible  phase  characterization  could  be  experi¬ 
mentally  obtained.  Recently,  a  theoretical  analysis  of 
the  open-circuit  capacitance  and  its  frequency  dependence 
has  been  reported  (Ref.  1) .  Measurements  in  7-mm  pre¬ 
cision  coaxial  line,  included  in  that  report,  were  con¬ 
sistent  with  the  theory,  but  not  quite  accurate  enough 
conclusively  to  discriminate  against  previously  reported 
values. 


The  HP  854 2 A  ANA  used  to  make  the  experimental 
measurements  for  Ref.  1  did  not  have  the  capability  of 
generating  synthesized  frequencies  referenced  to  a  quartz 
crystal  and  accurate  to  1  part  in  107  as  did  the  HP  8542B 
analyzer  used  for  this  research.  As  a  result,  a  more 
exact  experimental  corroboration  of  the  theoretical  values 
will  be  presented.  It  will  also  be  noted  that  the  rela¬ 
tive  frequency  variation  of  the  "effective  position"  is 
smaller  than  that  of  the  "effective  capacitance". 


Open  circuits  are  especially  convenient  high 
reflection  standards  when  used  along  with  short  circuits 
in  the  wide  band  (i.e.,  2-18  GHz)  calibration  of  ANA's 
since  the  attachment  of  several  different  offset  shorts 
(one  for  each  frequency  band,  2-4,  4-8,  8-12.4,  12.4-18 
GHz)  becomes  unnecessary:  the  result  being  a  significant 
reduction  in  wear  of  the  measurement  port  connectors. 
Additionally,  since  the  length  of  each  offset  short 
represents  the  A/4  plane  only  at  the  center  frequency  of 
each  band,  their  use  can  and  does  lead  to  noticeable 
band-edge  inaccuracies.  These  inaccuracies  are  illus¬ 
trated  in  Figs.  2-4,  which  represent  the  reflection 
measurement  of  an  APC-7  open  circuit  using  a  standard 
carefully  performed  GPM1  calibration  of  the  ANA.  Since 
it  is  reasoneib le  to  expect  the  reflection  of  an  open 
circuit  to  be  smooth  with  respect  to  frequency,  any 
discontinuities  would  have  to  result  from  the  calibration. 
As  can  be  seen,  when  the  marker  was  placed  on  each  of  these 
discontinuities,  the  frequencies  returned  by  the  analyzer 
were  indeed  4,  8,  and  12.4  GHz.  Use  of  the  open  circuit 
instead  of  offset  short  circuits  should  eliminate  these 
anomalies. 

Since  the  profits  to  be  gained  by  obtaining  an 
accurate  reflection  characterization  of  the  open  circuit 
seemed  so  substantial,  research  was  conducted  in  the 
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following  manner.  The  HP  8743A  Reflection  Test  Set 
(which  had  been  selected  for  extended  frequency  use  to 
18  GHz)  was  employed  throughout. 

The  trombone  delay  line  of  the  874 3A  was  adjusted 
so  that  the  time  delay  of  the  reference  path  was  sub¬ 
stantially  equal  to  that  of  the  signal  path.  This  was 
accomplished  by  manually  sweeping  the  8-12.4  GHz  frequency 
band,  and  adjusting  the  delay  line  until  the  display  of  the 
reflection  for  a  zero-plane  short  circuit  on  the  HP  8414 
Polar  display  approximated  a  dot.  This  precaution  mini¬ 
mizes  the  effect  of  frequency  jitter,  if  any. 

A  standard  HP  7-mm  calibration  kit  and  air  sliding 
load  provided  the  necessary  reflection  standards  'or  the 
calibration  of  the  ANA.  Previous  measurements  made  with 
this  calibration  kit  laid  suspicion  that  the  actual  con¬ 
structed  lengths  of  the  offset  shorts  differed  from  those 
specified  (Ref.  3) .  Therefore,  the  offset  shorts  were 
measured  using  a  Starrett  Model  65 3P  Dial  Comparator. 
Several  measurements  were  performed  around  the  center  con¬ 
ductor  of  each  offset  short  to  establish  that  the  shorting 
plane  was  indeed  perpendicular  to  the  line  (no  exceptions 
were  noted).  However,  Table  1  shows  that  the  nominal 
lengths  differed  from  the  measured  lengths  by  as  ch  as 
.12-mm.  As  a  result,  the  measured  lengths  were  used  in  all 
calculations. 


Offset 

Nominal 

Measured 

Difference 

Number 

Length  (mm) 

Length  (mm) 

(mm) 

4 

63.44 

63.32 

.12 

5 

31.70 

31.65 

.05 

6 

18.66 

18.58 

.08 

7 

13.36 

13.29 

.07 

Table  1 

-  Nominal  and 

measured  lengths 

of  HP  APC-7 

offset  shorts. 
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To  obviate  the  possibility  of  interaction  between 
zero-length  short  circuit  and  the  APC-7  bead,  offset 
shorts  were  used  for  both  the  reference  plane  and  the 
offset  reflection  plane.  Special  software  was  developed 
using  the  HP  error-correction  algorithms  but  which  allowed 
specification  of  the  offset  lengths  for  both  planes. 
Corrected  reflection  measurements  were  then  made  of  the 
open  circuit. 

Special  attention  was  given  to  frequencies  where 
the  difference  between  the  two  offset  lengths  was  nearly 
A/4,  where  calibration  accuracy  is  expected  to  be  optimum. 
For  each  of  these  frequencies,  four  sliding  load  measure¬ 
ments  were  taken,  each  with  the  load  displaced  by  A/8.  A 
complex  average  of  these  measurements  was  then  used  to 
establish  the  residual  reflection.  This  technique 
differed  from  HP's  which  uses  a  circle  fitting  routine  to 
find  the  residual.  In  some  cases,  where  the  load  ap¬ 
proaches  ideal  and  noise  is  present  in  the  measurements, 
HP's  circle-fitting  routine  could  diverge  and  thus  give  a 
defective  calibration.  The  complex  average,  however,  can 
never  yield  a  result  worse  than  the  largest  single 
measurement  in  the  average,  and  clearly  cannot  suffer  from 
inaccuracies  due  to  divergences.  Principal  weight  was 
given  to  each  of  the  favored  frequencies  mentioned  and  to 
the  fact  that  the  phase  of  the  open  circuit  must  be  zero 
at  zero  frequency  in  the  curve  fitting.  However,  a- 
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greement  at  other  frequencies  was  found  to  be  so  good  that 
these  precautions  may  not  have  been  necessary. 

Using  standard  linear  regression  techniques  to 
reduce  the  data,  the  following  frequency  characterization 
was  obtained,  where  A<J>  is  the  excess  (capacitive)  phase  of 
the  open  circuit,  in  radians,  and  f  is  frequency  in 
megahertz : 


A<)>  =  (5.02  x  10_5)f  +  (1.126  x  10"14)f3 


(1) 


This  result  implies  that  the  "effective  position"  of  the 
open  circuit  lies  beyond  the  physical  end  of  the  center 
conductor  by  a  distance  d(f)  in  millimeters  given  by: 


d(f )  =  1.198  +  2.69  x  10-10  f2 


(2) 


For  comparison  with  other  work,  the  effective 
capacitance  is  calculated  from 


2  7T  <fxlO°)Z. 


tan 


(3) 


where  ZQ  is  characteristic  impedance  of  the  transmission 
line  being  used. 
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The  empirical  results  for  the  effective  position  and 
effective  capacitance  are  tabulated  in  Table  2  along  with 
the  theoretical  results  given  in  Ref.  1,  where  the  interpo¬ 
lation  formula 


=  79.70 


-  ( f/34450) 


is  suggested. 

Table  2  shows  that  representation  by  effective 
position  (beyond  the  physical  end  of  the  center  conductor) 
is  less  frequency-dependent  than  the  effective  capacitance 
representation.  Thus,  the  effective  position  varies  by  7 
percent,  while  the  effective  capacitance  varies  18  percent, 
from  0-18  GHz. 

From  similitude,  the  effective  position  (or  capaci¬ 
tance)  for  14-mm  line  (at  half  the  frequency)  should  be 
twice  that  for  7-mm  line.  The  14  mm  GR900-WO  open  circuit 
has  a  closed  end,  however,  while  the  7  mm  open  circuits  in 
general  use  have  open  ends.  By  applying  similitude  to 
measurements  on  closed-end  7-mm  open  circuits,  we  find 
that  no  measurable  difference  in  effective  position  could 
be  attributed  to  this  constructional  difference.  Thus, 
the  effective  position  for  14-mm  line  would  range  from 
2.40  mm  at  1  GHz  to  2.57  mm  at  9  GHz.  This  range  falls 
within  the  specifications  for  the  GR900-WO  (2.40  -  2.80  mm) 
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and  leads  to  a  much  tighter  definition  of  the  14-mm  open 
circuit  them  one  could  hope  to  establish  with  a  precision 
slotted  line. 

The  empirical  characterization  of  the  open  circuit 
presented  here  was  undertaken  prior  to  the  publication  of 
the  theoretical  frequency  dependence.  Initial  close  a- 
greement  between  the  two  characterizations  was  so  im¬ 
pressive,  when  the  nominal  lengths  of  the  offset  shorts 
were  used,  that  we  were  inspired  to  measure  the  offset 
shorts.  The  substitution  of  these  measured  lengths  for 
the  nominal  lengths  reduced  the  difference  from  3  percent 
to  less  than  0.5  percent  between  the  two  characterizations. 

Also,  ANA  calibrations  ordinarily  include  a  flat 
short  circuit  placed  directly  at  the  connection  plane  of 
the  conductor.  The  characterization  of  the  open  circuit 
was  repeated  using  this  more  conventional  standard  along 
with  the  measured  offset  short  circuits.  Over  the  2-18  GHz 
range,  the  deviation  in  reflection  phase  between  this 
characterization  and  the  previous  one  was  nowhere  greater 
than  0.5  degrees  (See  Table  3).  This  observation  indi¬ 
cates  that  interaction  between  the  zero-plane  short  and 
the  measuring  port  bead  is  negligible. 

Based  on  these  results,  it  is  apparent  that 
frequency-corrected  open  circuits  can  replace  offset 
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Measurement  of 
open  circuit  phase 
based  on  zero-plane 


Frequency 

(MHz) 

short  circuit 

calibration 

(degrees) 

-A*, 

Eq.  1 
(degrees) 

2000 

-5.9 

-5.8 

3000 

-8.4 

-8.7 

4000 

-11.0 

-11.5 

6000 

-17.2 

-17.4 

8000 

-23.3 

-23.3 

10200 

-29.7 

-30.0 

12000 

-35.8 

-35.6 

12400 

-36.7 

-36.9 

14000 

-41.6 

-42.0 

14600 

-43.7 

-44.0 

16000 

-48.5 

-48.7 

18000 

-55.3 

-55.5 

Table  3  -  APC-7  open  circuit  phase  measurements  based  on 
zero-plane  short  circuit  calibration  compared 
to  empirically  derived  A<f>. 


shorts  as  high  reflection  standards  with  no  loss  of  phase 
accuracy,  compared  to  the  specified  accuracy  of  the 
HP  8542B  (reflection  phase,  from  1°  at  2  GHz  to  1.5°  at 


33 


ADAPTER  CORRECTION 

Since  the  HP  8542B  is  equipped  with  APC-7  precision 
connectors,  the  discussion  of  the  previous  section  shows 
how  a  wide-band  calibration  can  be  accomplished  without 
the  need  for  multiple  connections  of  offset  short  circuits 
in  this  primary  connector  system.  However,  there  are  so 
many  different  kinds  of  transmission-lines,  waveguides, 
and  connector  types  in  general  use  today  that  it  would  be 
entirely  impractical  to  construct  a  measurement  system 
based  on  each  type  available.  Therefore,  the  only  reason¬ 
able  solution  to  this  measurement  dilemma  is  to  perform 
these  measurements  through  passive,  reciprocal  adapters 
which  form  a  transition  between  the  primary  connector 
system  and  that  of  the  DUT  (secondary  connector  system) . 
The  present  method  for  carrying  out  measurements  in  the 
secondary  connector  system  requires  that  a  complete  cali¬ 
bration  be  made  (including  offset  shorts)  in  this  con¬ 
necting  system.  As  a  result,  this  technique  requires  that 
a  complete  calibration  kit  be  maintained  in  every  con¬ 
nector  and  transmission  format  used  in  measurements. 

The  cost  of  obtaining  and  maintaining  an  extensive 
inventory  of  these  calibration  kits  is  astounding  and  can 
be  prohibitive  when  appreciable  engineering  is  required  to 
design  and  prove  them.  For  example,  all  the  arguments 
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presented  in  the  previous  section  against  the  need  for 
connection  of  several  offset  shorts  in  APC-7  applies  even 
more  dramatically  to  SMA  lines  which  are  more  delicate  and 
easily  damaged.  Wear  of  these  connectors  is  an  important 
factor  and,  ideally,  the  adapters  and  standards  would 
have  to  be  replaced  at  frequent  intervals.  The  cost  of 
this  proposition  makes  it  impractical  and  unrealistic  for 
most  applications.  Aside  from  this  cost  factor,  any  pro¬ 
cedure  requiring  many  interconnections  for  each  calibration 
may  be  unreliable  when  wear  is  of  primary  concern. 

Of  course,  all  measurements  depend  upon  the  quality 
of  the  contacts  made  to  the  actual  DUT.  The  probability 
of  acquiring  an  accurate  calibration  would  definitely  be 
increased,  however,  if  fewer  contacts  are  made  during 
calibration.  The  technique  put  forward  in  the  previous 
section  of  using  an  open  circuit  characterization  to 
replace  the  offset  shorts,  is  equally  impractical,  because 
a  different  characterization  would  not  only  have  to  be 
made  for  each  connector  type,  but  also  for  each  sex  of 
type.  Pin  lengths  of  no  inherent  significance  to  the  con¬ 
nection  would  enter  into  some  of  the  open-circuit  charac¬ 
terizations  and  would,  therefore,  have  to  be  controlled. 

Further  complicating  the  problem  would  be  the 
questionable  propriety  of  mating  an  air  open-circuit 
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standard  to  a  dielectric-filled  connection  system.  A 
theoretical  characterization  would,  at  least,  be  more 
difficult  than  the  case  already  solved. 

For  these  reasons,  an  alternative  approach  of  ac¬ 
counting  for  the  adapter  calibration  was  adopted.  It 
employs  the  approximation  of  treating  the  adapter  as  dis¬ 
sipationless  (Ref.  2) .  This  approximation  is  not  at  all 
unreasonable  since  all  precision  adapters  are  constructed 
from  good  dielectrics  and  good  conducting  surfaces.  Under 
the  dissipationless  assumption,  one  needs  only  a  zero- 
reflection  standard  to  adjust  for  chart  center  and  a  zero- 
plane  short  circuit  to  define  phase  at  the  secondary  port 
connector . 

As  prescribed  by  this  adapter  correction  technique, 
the  ANA  is  calibrated  for  reflection  measurements  at  the 
primary  connector  port  according  to  the  methods  described 
in  the  previous  section.  Then  the  adapter  is  attached  and 
connected  with  what  is  assumed  to  be  an  ideal  termination. 
This  ideal  termination  can  be  simulated  by  a  precision 
fixed  load,  an  air  sliding  load  or  by  an  improved  tech¬ 
nique  where  a  line  and  fixed  load  are  computer  averaged  as 
discussed  in  a  later  section  under  Double  Running  Average. 
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For  each  frequency,  a  measurement  is  made  of  the 
thus-terminated  adapter,  corrected  with  respect  to  the 
primary  connector  calibration  and  stored  as  M^.  Then  the 
adapter  is  terminated  with  a  zero-plane  short  circuit,  the 
same  type  of  correction  measurement  is  accomplished  and 
stored  as  Mg.  Lastly,  reflection  measurements  are  made  on 
the  adapter  and  DUT,  likewise  corrected  with  respect  to 
the  primary  connector  port  and  stored  as  M^.  Then  the 
computer  performs  the  necessary  manipulations  to  correct  Mp 
for  reference  with  respect  to  the  secondary  connector  and 
represented  by  the  reflection  Fp. 


By  letting  S  represent  the  scattering  matrix  of  the 
adapter  where  port  1  is  the  primary  connector  and  port  2 
the  secondary,  then 


M_  -  S  +  S21  r° 

“d  “  sli  +  T-r-s 


F~ 
22  D 


(5) 


After  applying  the  principles  of  reciprocity  and  conser¬ 
vation  of  power  the  following  formula  for  adapter  cor¬ 
rection  is  derived  (Ref.  2) . 


r 


D 


MP  ~  “l  e“j29 
1  '  <  md  S 


(6) 
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where  m£  is  the  complex  conjugate  of  and  the  phase 

i  2  0 

factor,  e  J  ,  is  now  to  be  determined. 


With  the  short  circuit  as  a  phase  reference,  FD 

must  be  equal  to  -1.  Therefore,  substituting  M  for  M 

o  u 

and  -1  for  Fp  one  obtains: 


-1 


(7) 


or 


e 


-  j  2  9 


1  -  Ms  ml 

Ml  -  ms 


(8) 


where  Ml»  Mg,  Ml,  m£  are  all  known,  allowing  calculation 

i  2  0 

of  phase  factor  e  J 

If  El  represents  the  computer  calculation  of 


MS  Mt 


- 


(9) 


then  the  complete  complex  solution  for  the  adapter  cor¬ 
rection  becomes 


[El] 


(10) 
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Since  |e  J  |  must  always  equal  1,  and,  because 
limitations  on  significance  in  computer  calculations  may 
not  yield  a  result  where  ( El {  =  1,  El  is  normalized  by 
dividing  El  by  its  magnitude  | El |  before  applying  it  to 
the  final  adapter  correction. 

The  analysis  has  been  presented  for  a  single 
frequency.  The  computer  can  easily  perform  these  ma¬ 
nipulations  for  many  frequencies  in  little  time. 

Since  a  great  proportion  of  analyzer  measurements 
are  being  accomplished  through  adapters  to  other  con¬ 
nection  formats,  such  as  SMA,  this  technique  drastically 
simplifies  calibration  procedures  for  secondary  connector 
types.  In  addition,  it  reduces  the  chances  of  error  in 
calibrations  of  connector  types,  such  as  SMA,  which  are 
extremely  susceptible  to  lossy  connections,  by  requiring 
connection  of  fewer  standards  in  that  format.  The  ac¬ 
curacy  of  all  measurements  are  limited  by  the  quality  of 
the  DUT  connections,  but  calibration  errors  due  to  possi¬ 
ble  poor  connection  of  multiple  standards  is  minimized. 
Furthermore,  a  recalibration  is  not  necessary  when  the 
adapter  type  is  changed.  Only  the  measurement  of  the  new  a 
dapter terminated  by,  first,  its  matched  load,  and,  second, 
its  zero-plane  short  circuit  need  be  reaccomplished  in 
order  to  modify  the  adapter  corrections  for  the  new  con¬ 
nector  type.  This  feature  offers  the  advantage  of  a  great 
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time  savings  when  measurements  are  to  be  made  in  more  than 
one  connector  or  transmission  format  at  the  same  time. 

Considerable  interchanging  of  adapters  is  necessary 
for  the  measurement  of  non-insertable  two-ports  (e.g., 
devices  with  two  connectors  of  the  same  sex) .  The  ap¬ 
plication  of  the  adapter  correction  technique  to  this  case 
will  be  discussed  under  Extension  of  Adapter  Correction  to 
Two-Port  Device  De-Embedding. 
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CHARACTERIZATION  OF  PERFECTLY  MATCHED  TERMINATION 

In  the  absence  of  a  perfectly  matched  termination 
in  the  real  world,  beadless  air-dielectric  sliding  loads 
have  been  the  standard  of  choice  for  calibrating  the 
HP  8542B  ANA.  Three  or  more  measurements  must  be  made  on 
the  sliding  load  at  each  frequency.  Then  a  circle  fitting 
algorithm  (CENT)  must  be  applied  to  find  the  center  of  the 
complex  circle  circumscribed  by  these  points  (the  center 
representing  the  reflection  of  an  ideal  load) (Ref.  6) . 

There  are  at  least  four  major  drawbacks  to  this  method  of 
characterization. 

First,  in  using  a  sliding  load  for  measurements 
over  wide  frequency  bands,  three  positions  will  not  suf¬ 
fice  for  accuracy.  Four  or  five  measurements  at  each 
frequency  are  required  and  choosing  the  length  of  each 
slide  so  that  the  measurements  form  a  relatively  well 
defined  circle  at  all  frequencies  is  no  trivial  matter. 
Second,  the  CENT  program  used  for  the  circle  fitting  can 
diverge  and  introduce  large  errors,  when  used  to  fit  data 
collected  on  a  load  which  approaches  perfect  at  some 
frequencies  and  which  may  be  displaced  by  noise  in  the 
measurements.  The  above  problems  have  been  studied  (Ref.  7) . 
That  is,  the  optimum  pattern  of  slides  to  minimize  the 
errors  introduced  by  the  CENT  program  has  been  sought  for 


broad-band  calibration.  In  addition,  of  course,  one  must 
expect  errors  in  the  measurement  process. 

Third,  this  established  method  of  calibration  makes 
no  use  of  information  at  other  frequencies  to  assist  in 
the  calibration  at  any  particular  frequency.  Fourth, 
some  connection  systems  (e.g.,  the  popular  SMA  system) 
are  based  on  the  concept  of  complete  filling  with  a  solid 
dielectric  (e  %  2.08)  throughout  the  connectors  and 
transmission  line.  The  use  of  an  air-dielectric  sliding 
load  as  a  standard  in  a  solid-dielectric  transmission 
format  can  be  expected  to  lead  to  measurable  error  at 
relatively  high  frequencies  within  the  operating  range  of 
this  connection  format.  Ideally,  the  reflection  coef¬ 
ficient  of  a  standard  should  be  specified  with  respect  to 
the  transmission  system  used  to  interconnect  the  com¬ 
ponents  or  DUT. 

The  last  problem  could  be  addressed  by  applying 
a  correction  when  an  air-dielectric  standard  is  used.  The 
diameter  changes  implied  by  a  dielectric  constant  can  be 
expected  to  lead  to  a  discontinuity  capacitance  at  the 
transition  from  dielectric  to  air.  This  discontinuity 
capacitance  can  be  theoretically  estimated  and  applied  as 
a  small  correction,  its  importance  increasing  with 
frequency.  Our  experimental  investigations  of  this  matter 
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quickly  showed  that  practical  connectors  introduced  much 
larger  discrepancies  (with  more  complicated  frequency 
dependence)  than  the  estimated  discontinuity  capacitance. 
As  a  result,  it  was  considered  all  the  more  important  to 
be  able  to  use  practical  terminations  as  low  reflection 
standards,  when  evaluating  devices  intended  for  use  with 
such  interconnecting  systems. 

The  policy  of  using  a  physical  transmission  line 
and  a  fixed  load  as  a  standard  of  zero  reflection  has  long 
been  advocated  for  swept- frequency  measurement  of  re¬ 
flection  magnitude  (Refs.  8,  9,  10).  These  methods  only 
had  the  capability  of  dealing  with  reflection  magnitudes 
and  were  not  capable  of  extracting  phase  information  from 
the  measurements.  On  the  contrary,  implementation  of  a 
similar  policy  on  the  ANA  is  not  only  convenient,  but  is 
also  capable  of  performing  and  recording  phase  measure¬ 
ments  . 


A  time-domain  ref lectometer  has  been  synthesized 
in  software  (TIMED)  by  Fourier  transformation  of  re¬ 
flection  data  obtained  in  the  frequency  domain  with  an 
ANA  (Ref.  11) ,  In  the  time  domain,  the  reflection  from 
an  imperfect  termination  of  a  finite  length  of  transmission 
line  can  easily  be  identified.  One  could  replace  this 
time-domain  reflection  by  an  extrapolation  equivalent  to  a 


perfect  termination.  The  inverse  Fourier  transform  would 
then  provide,  at  all  frequencies,  the  equivalent  of  cali¬ 
bration  measurements  on  a  perfect  load. 

The  discussion  of  TIMED  indicates  that  wide-band 
measurements  of  a  fixed  length  of  imperfectly  terminated 
transmission  line  should  contain  the  information  necessary 
to  characterize  a  perfectly  matched  termination.  However, 
the  present  research  has  established  the  validity  of  a 
simpler  method,  using  only  frequencies  in  the  general 
vicinity  of  the  particular  frequency  of  interest  (frequen¬ 
cy  window) . 

Recently,  the  development  of  3.5  mm  air-dielectric 
line  (WSMA)  was  purported  to  relieve  "nagging  SMA 
measurement  problems"  (Ref.  12) .  However,  this  approach 
still  does  not  avoid  the  compromise  of  defining  a  zero 
reflection  standard  for  a  dielectric  environment  using  air 
lines  as  standards.  The  method  to  be  presented  here 
provides  a  practical  solution  to  all  of  the  four  problems 
mentioned  above. 

The  errors  in  the  reflection  measurement  are 
assumed  to  be  linear,  as  in  the  present  calibration  pro¬ 
cedures  for  the  HP  8542B.  The  ultimate  goal  is  to  determine 
the  uncorrected  reading  that  would  be  given  by  the  system 
if  a  perfect  load  was  obtainable:  "perfect"  here  meaning 
a  match  to  the  intended  connecting  system.  Since  no 


"perfect  load"  is  physically  available,  this  is  ac¬ 
complished  by  making  a  corrected  measurement  of  reflection 
for  a  practical  load.  This  reading  will  be  referred  to  as 
the  reflection  residual  of  the  load  or,  simply,  residua 
(Tr).  After  these  measurements  are  acquired  for  each 
frequency,  they  can  be  applied  as  corrections  to  measure¬ 
ments  on  DUT's  thereby  synthesizing  a  "perfect  load" 
reference. 

The  method  of  characterizing  a  perfect  load  (without 
resorting  to  Fourier  analysis)  is  based  on  the  assumption 
that  the  system  errors  vary  relatively  slowly  with  frequen¬ 
cy.  Relatively  slowly  means  that  the  round-trip  re¬ 
flection  delay  time  T  for  the  transmission  line  being  used 
as  a  standard  is  such  that  the  variation  in  the  system 
error  may  be  neglected  over  a  frequency  interval  1/T. 

Also  assumed  is  that  the  imperfect  termination  has  a  re¬ 
flection  coefficient  whose  magnitude  also  varies  rela¬ 
tively  slowly  with  frequency  in  the  same  sense.  The  re¬ 
flection  coefficient  of  the  load  is  considered  to  include 
that  of  the  connector  used  to  attach  it  to  the  end  of  the 
transmission  line.  Even  with  this  inclusion,  the  frequen¬ 
cy  dependence  of  the  physical  load  would  be  expected,  on 
the  grounds  of  small  size  and  nonresonant  design,  to  be 
less  of  a  problem  than  the  system  errors. 
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The  input  reflection  coefficient  of  such  a  termi¬ 
nated  line  of  electrical  length  l  would  be 


(11) 


where  a  is  the  attenuation  constant,  6  is  the  propagation 
constant  (R  *  f  where  R  =  2.0965  x  10  4  radians/MHz/cm  in 
air  and  f  is  frequency  in  MHz) »  and  is  the  reflection 
coefficient  of  the  imperfect  load.  Substituting  R  *  f 
for  0  and  rearranging  Eq.  (11)  yields 

r  *  e“2a£  e"2jR£f  r.  (12) 

Xj 

If  the  frequency  dependence  of  a,  R,  and  rL  can  be  neg¬ 
lected,  the  input  reflection  will  describe  a  circle  as 
frequency  is  varied.  The  average  of  this  complex  re¬ 
flection  over  a  frequency  interval  1/T  (frequency 
window) ,  being  the  center  of  this  circle,  will  vanish. 

However,  because  the  ref lectometer  used  for 
measurement  is  not  perfect,  the  input  reflection  will 
still  trace  a  circle,  but  the  center  of  the  circle  will 
be  offset  from  zero  due  to  system  errors.  The  input  re¬ 
flection  can  be  represented  as 
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r  =  rR  +  a  e"2jR£f  rL  (13) 

j  6 

where  A  =  e_2a£  and  T.  =  I T_  |e  L.  If  the  frequency 

dependence  of  TL  is  neglected  over  one  cycle,  @L  is  a 

constant,  and  T  is  the  center  of  the  circle  representing 
R 

the  system  error  and  is  usually  accepted  as  being  an 
approximation  to  the  uncorrected  measurement  of  a  perfect 
load. 


The  real  and  imaginary  components  of  f  become 


Re  [ T ]  =  Re  [TR]  +  A| rL f  cos  (2Rif  -  0L>  (14) 

and 

Imag  [f]  =  Imag  (fR]  -  A  |  |  sin  (2RJtf  -  0L>  (15) 


which  are  clearly  the  sum  of  a  constant  (reflection 
residual  or  system  error)  and  a  sinusoidal  component 
(ripple  factor)  caused  by  the  imperfect  load. 

Since  the  network  analyzer  is  capable  of  making 
complex  reflection  measurements,  a  rather  simple  algo- 
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rithm  can  be  used  to  average  out  the  sinusoidal  com¬ 
ponents.  This  running  average  is,  therefore,  equivalent 

★ 

to  applying  a  low-pass  filter  yielding  an  adjusted 

r  =  r_  which  can  now  be  used  to  correct  for  the  re- 
a  R 

flection  residuals  (system  error) . 

Two  variations  of  this  running  average  technique 
were  investigated  for  this  project  and  each  will  be 
described  in  detail  shortly.  A  Single  Running  Average 
Technique  was  effective,  but  required  a  minimum  of  23 
points  per  cycle  for  the  average.  However,  a  Double 
Running  Average  became  the  filter  of  choice  for  the  new 
HP  8542B  calibration  technique,  because  it  offered  the 
same  effectiveness  as  the  Single  Running  Average  while 
only  requiring  10  points  to  compute  the  average. 

The  initial  attempts  to  apply  the  above  procedures 
(calculating  the  Single  Running  Averages  of  the  real  and 
imaginary  parts  of  the  uncorrected  reflection  coefficient) 
did  not  work  well  at  all.  The  system  errors,  particularly 
in  band,  did  not  satisfy  the  condition  of  being  slowly 
varying  over  the  168  MHz  frequency  interval  required  by 
the  reference  line  of  89.4  cm.  electrical  length  used. 

The  prospect  of  using  a  longer  line  and  closer  spaced 
frequencies  was  unattractive. 


low-pass  in  the  time  domain 
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Fortunately,  it  was  conjectured  subsequently  that 
these  variations  must  be  caused  by  echos  within  the  long 
lines  inside  the  analyzer  and  were,  therefore,  a  part  of 
every  uncorrected  or  raw  measurement.  Consequently, 
complex  division  of  any  raw  measurement  by  the  raw  re¬ 
flection  measurement  of  a  short  circuit  should  eliminate 
these  variations  greatly  reducing  the  frequency  dependence 
of  the  errors.  In  fact,  this  relatively  simple  partial 
correction  of  normalizing  by  division  did  suppress  the 
apparent  echoes  and  revealed  the  sinusoidal  periodicity 
well  enough  to  conduct  a  plausible  running  average.  One 
can  remultiply  the  averaged  quotient  by  the  raw  measure¬ 
ment  of  the  short  circuit  to  establish  the  "raw  measure¬ 
ment"  one  would  have  obtained  if  the  reference  trans¬ 
mission  line  was  terminated  perfectly.  Of  course,  any 
high-reflection  standard  could  be  used  for  the  normal¬ 
izing.  This  concept  afforded  the  possibility  of  cali¬ 
brating  the  ANA  without  the  need  for  a  sliding  load. 

Figure  5  shows  the  real  part  of  a  normalized  re¬ 
flection  coefficient  measurement  of  a  terminated  SMA  line 
with  em  electrical  length  of  89.4  cm.  superimposed  by 
the  single  running  average.  Normalization  was  accomplished 
by  dividing  the  uncorrected  measurements  by  the  uncor¬ 
rected  measurements  of  a  short  circuit  placed  at  the  SMA 
end  of  the  adapter.  The  Single  Running  Average  is  based 
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Fig.  5.  Results  of  Single  Running 
Average  using  23  points 
per  cycle. 


on  23  points  equally  spaced  in  frequency.  Thus,  the 
imperfections  of  the  load  are  manifested  in  the  variations 
having  a  period  of  168  MHz  in  the  frequency  domain. 

It  was  initially  hoped  that  this  technique  of 
employing  a  running  average  could  be  adopted  as  a  general 
replacement  for  sliding  loads,  but  subsequent  experi¬ 
mentation  revealed  that  for  APC-7  lines  with  their  excel¬ 
lent  impedance  characteristics,  precision  connectors,  and 
air  environment,  the  established  method  of  employing  a 
sliding  load  for  characterization  could  not  be  improved 
upon.  Therefore,  the  sliding  load  technique  was  retained 
for  APC-7  calibrations.  However,  for  connection  systems 
such  as  SMA,  where  the  dielectric  filling  and  the  extreme 
fragility  of  sliding  loads  makes  the  whole  concept  of 
referencing  against  transmission  lines  more  relevant,  the 
results  were  remarkably  satisfying.  As  a  result,  the 
averaging  technique  was  adopted  to  characterize  a  perfect 
load  for  the  adapter  correction  sequence. 

Hence,  the  following  calibration  procedures  for  the 
HP  8542B  were  developed.  The  system  is  first  calibrated 
in  APC-7  using  the  short,  open,  and  sliding  load  as 
standards  at  all  measurement  frequencies  for  APC-7  or  at 
all  frequencies  necessary  to  oerform  a  running  average  for 
adapter  corrected  measurements.  If  adapter  correction  is 


called  for,  then  measurements  are  made  on  the  adapter  and 
short  and  then  on  the  adapter  and  transmission  line 
reference.  These  measurements  are  then  corrected  against 
the  APC-7  calibration.  Since  this  step  corrects  the  trans¬ 
mission  line  measurements  for  system  errors,  it  became 
unnecessary  to  perform  the  normalization  described  earlier 
while  still  attaining  the  same  results.  The  Double 
Running  Average  is  then  computed  and  stored. 

Measurements  are  then  carried  out  on  the  adapter 
and  DUT  at  the  ultimate  frequencies  of  interest,  corrected 
against  the  APC-7  calibration,  and  then  adapter  corrected. 
The  results  are  output  in  any  of  a  variety  of  formats. 
Repeat  measurements  can  be  readily  performed  on  the  same 
adapter  without  recalibration. 

Another  important  factor  is  that  when  employing 
a  computer  with  sufficiently  large  storage  capacity,  a 
change  of  adapter  can  be  achieved  by  simply  redoing  the 
adapter  correction  sequence.  Present  procedures  for 
HP  8542B  require  a  complete  calibration  in  the  new  trans¬ 
mission  format,  an  extremely  tedious  operation  with  ever 
present  hazards  of  normal  wear,  or  of  abnormal  wear  on 
adapter  or  standards  from  misconnection  or  improper 
tiv  tening. 
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The  new  calibration  procedures  described  here  offer 
the  following  advantages.  The  number  of  standards  in  the 
APC-7  cal  kit  are  reduced,  calibration  kits  in  other  trans¬ 
mission  formats  become  unnecessary,  sliding  loads  other 
than  APC-7  are  eliminated,  and  operator  interaction  is 
reduced.  The  benefits  far  outweigh  the  relatively  minor 
drawbacks  of  slightly  longer  computation  time  and  the 
increased  computer  memory  (a  very  inexpensive  commodity 
today) . 
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SINGLE  RUNNING  AVERAGE 

The  Single  Running  Average  was  the  first  filtering 
technique  investigated  for  this  study.  Even  though  it  was 
eventually  abandoned  in  favor  of  the  Double  Running  Aver¬ 
age  (to  be  discussed  later) ,  valuable  information  was 
gained  during  the  experimentation  which  was  directly  bene¬ 
ficial  to  the  application  of  the  Double  Running  Average. 

It  is  of  significance  to  note  here  that  the  Single  Running 
Average  was  not  replaced  because  it  did  not  function  well, 
but  because  it  required  a  minimum  sample  size  of  23  points 
per  cycle  in  order  to  work  adequately. 

Figure  5  shows  the  corrected  real  component  of  a 
reflection  measurement  of  a  terminated  89.4  cm  SMA  line 
superimposed  with  its  Single  Running  Average.  Twenty- 
three  samples  per  cycle  were  used,  and,  as  can  be  seen, 
the  ripple  factors  are  diminished  to  an  extent  that  they 
are  rendered  invisible.  In  Fig.  6,  the  real  component  for 
the  same  line  is  presented  superimposed  by  a  Single 
Running  Average  with  a  sample  size  of  21  points  per  cycle. 
The  ripple  factor  has  clearly  become  visible  in  this  plot. 
After  thorough  investigation,  the  conclusion  was  reached 
that  a  minimum  of  23  samples  per  cycle  were  required  for 
the  single  running  average. 
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Results  of  Single  Running 
Average  using  21  points 
per  cycle. 
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The  Single  Running  Average  was  set  up  and  computed 
in  the  following  manner.  A  frequency  window,  AF,  was 
chosen  such  that  2RH  AF  =  2  it  and  then  the  input  reflection  F 
is  measured  at  n  equispaced  points,  f^,...,f  ,  separated 
by  a  frequency  interval,  Af,  where  n  ^  23  and  odd  such 
that  (n-l)Af  =  AF.  The  real  and  imaginary  components  of  T 
are  then  averaged  over  AF  weighting  the  redundant  end 
points  (Fig.  7a).  These  end  points  were  weighted  by  1/2, 
and  their  combined  total  treated  as  one  point  to  obviate 
the  chance  of  biasing  caused  by  the  redundancy.  The  si¬ 
nusoidal  ripple  factors  vanish  leaving  an  adjusted  T  =  TR 
for  f  the  center  frequency  of  AF  (Fig.  7b)  . 

Two  separate  and  distinct  methods  of  employing  the 
Single  Running  Average  were  developed,  each  having  its  own 
advantages  under  different  circumstances.  The  first 
method,  called  discrete  running  average,  was  set  up  by 
simply  measuring  F  at  --^-Af  frequency  intervals  each  side 
of  the  frequencies  of  interest  (FI's)  and  applying  the 
algorithm  exactly  as  described  above.  This  method  was 
particularly  suited  for  wide  band  measurements  where  the 
FI's  are  far  apart  (where  difference  between  adjacent 
FI's  >  AF) .  Indeed,  the  discrete  running  average  was 
adaptable  to  most  any  interval  of  FI's,  but  suffered  from 
the  following  disadvantages.  First,  a  minimum  of  23  dis- 


Plot  of  Re[r]  (a)  before  and  (b)  after 
averaging. 
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tinct  measurements  were  required  for  each  FI  resulting  in 
the  need  for  large  amounts  of  computer  storage  for  compu¬ 
tation.  Also,  since  the  complete  algoithm  had  to  be  com¬ 
puted  for  each  frequency  of  interest,  the  algorithm  was 
rather  slow. 

A  second  method,  called  a  sliding  running  average, 
was  conceived  which  could  handle  the  special  case  where  the 
difference  between  adjacent  FI's  =  A f.  Thi'  algorithm 
was  particularly  sjiiable  for  relatively  narrow-band 
measurements  and  offered  the  advantages  of  reduced  storage 
requirements  and  remarkably  faster  speed. 

For  the  first  FI  (FI^) ,  the  discrete  method  was 
used  to  compute  the  running  average.  However,  since  the 
difference  between  adjacent  FI's  equals  Af,  the  calcu¬ 
lation  of  the  algorithm  for  FI_,...,  FI  was  much  simpler. 

c  m 

Since 


rf  /2  +  rf  +  rf  +....+  rf  +  rf  /2 

1 _ r2  3 _  n-l  n 

n  -  1 


(16) 


and 


/2  +  r-  +  +  ....  r-  +  r,  /2 

_2 _ r3  r4  n  £n+l 


n-l 


(17) 


Then  T 


__  could  be  quickly  calculated  by  subtracting  =-7  -?r 

FI2  rf.  rf  rf  _  2(n"1) 


and 


n+1 


JTn^rr  and  bY  addin9  rn^ir  and  ruPir  to  rpi1- 


For  FI ^  , « •  • ,  FIm,  the  process  is  repeated  in  the  same 
manner.  This  process  eliminated  the  need  for  storing 
measurements  at  overlapping  frequencies  and  greatly  reduced 
the  number  of  calculations  required  for  second  and  suc¬ 


ceeding  rFI's. 


The  disadvantages  of  this  method  were  that  it  was 
only  adaptable  to  the  special  cases  where  measurements 
are  being  made  over  a  relatively  narrow  band  and  that  the 
frequency  spacing  was  limited  to  Af. 


Since  the  Single  Running  Average  required  such  a 
large  sample  size  to  compute,  investigations  into  another 
running  average  technique  ensued.  Several  aspects  of  the 
discrete  and  sliding  average  proved  to  be  useful  in  the 
development  of  the  Double  Running  Average  Technique,  to  be 
discussed  next. 
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DOUBLE  RUNNING  AVERAGE 

Because  the  number  of  points  per  cycle  necessary  to 
compute  the  Single  Running  Average  was  so  large,  it  was 
hoped  that  another  digital  low-pass  filter  could  be  found 
which  was  much  more  efficient  (i.e.  required  fewer  points 
per  cycle) . 

A  computational ly-simple  candidate  was  the  Double 
Running  Average.  First,  the  Single  Running  Average  is 
applied  to  the  corrected  real  and  imaginary  measurements 
on  the  load.  Then  the  resulting  data  is  Single  Running 
Averaged  again  to  obtain  the  final  result  at  the  frequency 
of  interest.  Thus,  the  corrected  real  and  imaginary 
measurements  are  subject  to  a  Double  Running  Average.  This 
process  is  equivalent  to  a  single  weighted  average  with  a 
triangular  set  of  weights,  but  avoids  multiplications. 

This  procedure  also  had  the  gratifying  effect  of 
reducing  the  number  of  measurement  points  required  per 
cycle  of  ripple  from  23  to  5  for  an  equally  satisfactory 
suppression  of  the  ripple.  However,  as  a  factor  of 
safety,  7  measurements  per  cycle  were  adopted  for  the 
technique. 


60 


Mathematically,  the  Double  Running  Average  is 
computed  in  the  following  manner.  Treatment  of  the  real 


component  will  be  described.  The  imaginary  component 
is  found  in  exactly  the  same  way. 


A  frequency  window  of  2AF  must  be  used  where 
2R£  AF  =  2rr  (frequency  window  must  be  two  cycles  wide)  . 
The  system-corrected  measurements  of  the  reflection  A (I) 
are  obtained  at  4M  +  1  equispaced  points,  f^,..., 
f2M+l'***'  f4M+l'  seParated  by  a  frequency  interval,  Af, 
where  M  >_  3  such  that  2MAf  =  AF,  and  where  ^2m+1 
FI.  Let  B ( I )  =  Re [A( I) ]  at  each  frequency,  then  for 
J  =  M  +  1  to  3M  +  1  let 


D(J,  =  +  3T  B<K)  +  fiijisa. 

*  K=J-M+1  1 


As  in  the  case  for  the  Single  Running  Average  the  re¬ 
dundant  end  points  are  weighted  by  1/2  and  their  combined 
weight  is  treated  as  one  point  in  the  average. 


The  second  average  is  now  computed  yielding 


2<S=SL  ♦  "T1  dip)  ♦  EiSiSl- 


P=N-M+1 

36 


.  ,  '“"i— i.  '-Jvm* 
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where  N  =  2M  +  1  or  f.,  is  the  FI.  As  can  be  seen,  the 
redundant  end  points  are  weighted  and  treated  in  the  same 
manner  as  above.  Figure  8  shows  a  pictorial  view  of  the 
progression  of  the  averages. 

As  in  the  case  of  the  Single  Running  Average,  there 
are  two  different  methods  of  employing  the  Double  Running 
Average,  discrete  and  sliding.  Of  particular  significance, 
however,  is  that  in  both  the  discrete  and  sliding  Double 
Running  Average,  the  first  application  of  the  Single 
Running  Average  lends  itself  well  to  the  use  of  the  Sliding 
Running  Average.  A  considerable  reduction  in  memory 
fetches  and  computation  time  are  realized  by  taking  ad¬ 
vantage  of  this  simplification. 

In  the  case  where  measurements  are  being  made  in 
such  a  way  that  ultimate  FI's  in  the  Double  Running 
Average  are  equispaced  Af  apart,  the  sliding  average 
principal  can  be  applied  to  both  the  first  and  second 
Single  Running  Averages  resulting  in  an  even  greater 
savings  of  computation  time  and  an  extra  benefit  of  re¬ 
duced  computer  memory  requirements  compared  to  the  discrete 
Double  Running  Average  technique.  This  technique,  however, 
is  subject  to  the  same  limitations  discussed  under  the 
sliding  Single  Running  average  of  being  only  useful  over 
relatively  narrow  bands  and  of  having  a  frequency  spacing 
limited  to  Af. 


Diagram  showing  progression  of  averages  in 
Double  Running  Average. 
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. .  <■ 


Adoption  of  the  Double  Running  Average  over,  r >e 
Single  Running  Average  resulted  in  a  70%  reduction  in  the 
required  number  of  points  per  cycle  while  only  increasing 
the  number  of  computations  by  approximately  23%.  On  a 
computer  with  only  8K  of  memory,  the  decrease  in  storage 
requirements  was  a  much  more  important  factor  than  compu¬ 
tation  time.  In  fact,  investigations  showed  that  the 
computation  time  was  fast  enough  that  the  increase  was 
hardly  noticeable. 

The  minimum  number  of  samples  per  cycle  required 
for  the  Double  Running  Average  was  experimentally  found 
to  be  five  (Fig.  9) .  This  number  is  only  slightly 
larger  than  the  theoretical  limit  of  three  imposed  by  the 
Sampling  Theorem  for  a  periodic  waveform.  In  most  practi¬ 
cal  applications  of  this  theorem,  the  limit  is  usually 
accepted  to  be  twice  the  theoretical  limit  as  a  factor  of 
safety.  Therefore,  the  seven  samples  per  cycle  used  for 
the  Double  Running  Average  appeared  to  be  the  smallest  odd 
number  of  samples  which  could  reasonably  be  expected  to 
perform  successfully  in  a  practical  sense. 

Since  the  Double  Running  Average  did  an  equally  good 
job  of  removing  the  ripple  factor  as  the  Single  Running 
Average  did  (see  Figs.  9  and  10)  and  since  further  re¬ 
duction  in  the  number  of  sample  points  per  cycle  seemed 
unpromising  on  the  basis  of  the  above  reasoning,  it  became 
the  digital  filter  of  choice  for  this  project. 
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Fig.  9.  Results  of  Double  Running 
Average  using  5  points 
per  cycle. 
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DERIVATION  OF  CALIBRATION 

For  all  the  reflection  calibrations,  the  standard 
HP  linear  error  correction  model  (Ref.  13)  was  used  which 
is  diagrammed  below: 


Fig.  11.  HP  error  correction  model  for 
reflection  measurements. 

and  which  assumes  that  *  S2^  =  0.  Using  this 

REFL 

flowgraph  to  solve  for  ==■ =—  =  r  the  following  relation 

rvbr 

is  extracted 


r  =  e  + 

MEAS  00  1-e 


11  S11 


(20) 


where  rMEAS  represents  the  measured  reflection  coefficient 
of  DUT  and  S^,  the  actual  reflection  coefficient  of  DUT. 
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eQ0  is  the  crosstalk  error  due  to  imperfect  directivity 
of  the  coupler,  e^  is  the  error  due  to  the  imperfect 
gain  tracking  between  test  and  reference  channels,  and  e^ 
is  the  mismatch  error  due  to  imperfect  Port  1  match. 

It  is  clear  that  three  reflection  standards  could 
be  used  to  solve  for  the  three  error  correction  coef¬ 
ficients  eQ0,  eQ1,  and  e^.  The  three  standards  chosen 
for  this  project  were  the  short,  open,  and  character¬ 
ization  of  a  perfectly  matched  load.  If  Tt  is  the 
measured  reflection  coefficient  of  the  load,  rg  of  the 
short,  and  of  the  open  circuit,  then  S^fL)  =  0, 

S1;L(S)  =  -1,  and  Sn(0)  =  V  where  r  =  le“jt,)(f).  It  has 
been  shown  that  for  the  7-mm  precision  50  ohm  open-circuit 
<J>(f)  =  5.02  x  10  ^  f  +  1.126  x  10  f3  radians  where  f  is 
frequency  in  MHz.  Substituting  the  stated  values  of 
into  Eq.  (20)  for  each  standard,  one  obtains  the  following 
equations : 

rL  =  e00  (21) 

r  =  e - fO-1- 

S  e00  1  +  en 


(22) 
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r 


o 


'00 


(23) 


Substituting  TL  for  eQg  in  (22)  and  (23)  and  rearranging 
terms  gives 


r  -  r  =  "  6qi 

S  L  1  4  e 


11 


r  r  -  e°l  F 
0  L  '  1  -  eu  t 


(24) 


(25) 


Dividing  Eq.  24  by  25  and  multiplying  both  sides  by  T,  one 
finds 


r 


1 

I 


(26) 


Let 


Q  = 


(27) 


which  is  now  calculable  and  Eq.  26  becomes 


Q 


ell  F  ~ 
ell  + 


1 


(28) 


From  Eq.  22 

601  *  (PL  *  rs)(1  +  ell>  (30) 

The  system  now  has  enough  information  to  make  fully 
corrected  reflection  measurements  (is  "calibrated") . 
Equation  20  can  now  be  reversed  to  find  in  terms  of 

^MEAS  wll^ch  is  the  ultimate  aim  and  it  becomes 


Sn  (DUT) 


(31) 


Therefore,  fully  corrected  reflection  measurements  can  be 
made  after  first  calibrating  the  system  using  the  three 
reflection  standards  with  known  S-^'s. 


CALIBRATION  AND  MEASUREMENT  PROGRAMS 


The  calibration  and  measurement  programs  will  now 
be  presented  in  detail.  The  programs  are  written  in  HP 
combined  ANA,  ATS,  and  TODS-II  BASIC  languages,  and  run  on 
an  HP-2100S  computer  which  is  equipped  with  8K  usable 
memory,  and  convenient  graphics  output  devices. 

Due  to  core  memory  limitations ,  the  complete 
sequence  (including  adapter  correction)  is  controlled  by 
six  programs  which  are  maintained  in  disc  files  and  auto¬ 
matically  chained  into  memory.  The  APC-7  sliding  load, 
short,  and  open  are  used  for  the  initial  calibration.  The 
adapted  zero-plane  short  and  running  averaged  load  are 
employed  in  the  adapter  correction.  The  discrete  Double 
Running  Average  is  incorporated  for  maximum  flexibility. 

Figure  12  shows  the  flowgraph  for  program  1  which 
is  tasked  with  making  all  measurements  necessary  for  the 
APC-7  calibration.  Figures  13-14  show  the  listing  of 
program  1.  It  first  requests  the  electrical  length  of  the 
reference  line  which  will  be  used  for  the  load  character¬ 
ization  in  the  adapter  correction.  From  this,  it  computes 
F  =  AF  and  displays  the  minimum  and  maximum  frequencies 
available  within  equipment  limits,  and  asks  for  the  start, 
stop,  and  step  frequencies  in  order  to  compute  the  FI's. 

It  then  computes  F9  =  Af  and  N(the  number  of  FI's),  F3 
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PfiSE  1 


1  PEH  THIS  PROGRAM,  DEVELOPED  ON  IS  DESIGNED 

2  PEN  Tu  MAKE  AN  APC-7  CALIBRATION  IN  PREPARATION  FOP 

2  PEM  THE  DOUBLE  RUNNING  AVERAGE  IN  THE  ADAPTEP  COPPECTION 
4  PEM  REV  1  1  -  2  4 ' 8  0  REM  1  2  f  2  2  ✓  3  0  1  :  A  N  U  7  2  .  S 

2  0  COM  2  C 1 3 , 511. LC13, 513, 0113. 313 
40  COM  F  C  5  1  J  ,  N  t  1  0  J 

1O0  D  C  P  L  A  Y  “ U  H  A  T  IS  L  E  N  G  T  H ( C  M )  OF  L I N  E " j 
120  BELL 
140  INPUT  L 

130  PEM  CALC  1  -PER IOD  FOR  LENGTH  OF  LINE  AND  UPPER  AND 

131  PEM  LOWER  FREQ  LIMITS  FOR  RUN  AVG 
13  0  LET  F«2?-?S0.(2*L) 

ISO  LET  F  1  »  I  N  T  t  101+F) 

2O0  L  E  T  F  2 • I N T  (  1 3  0  0  0 - F l 

2  20  D  2  P  L  A  Y  “FREQ  (  MHZ!  AVAILABLE  -  M  I  N  .  “  ;  F  1  ,  "  M  A  X  •  “  ;  F  2 
24  0  D I PLAV  “FREQ  <  MH2 )  -  START. STOP, STEP “( 

230  BELL 

2  3  0  input  f  i  ,  f  2  ,  s 

230  PEM  CALC  NUMBER  OF  FREQ  STEPS  ASKED  FOR  IN  THE  OUTPUT 
200  LET  N.lNTH*iF2-Fl)/S) 

3lO  PEM  INITIALISE  THE  MEASUREMENT  EQUIPMENT 

320  LET  F  9  «  F  ^  3 

240  LET  F3.F1-G4F9 

3  3  0  F  C  A  L  f ( F  3  I 
380  BCNT  1  ( F3  i 

4  0  0  :  7  E  L  1  I  1  1  : 

420  WAIT  i30i 

430  PEM  CALC  AND  STORE  THE  OUTPUT  FREQS 

440  FOR  K « 1  TO  N 

480  LET  FIKI.F1* (K-l ) *S 

430  NEXT  K 

300  DSPLAY  “CONNECT  APC-7  SHORT” 

320  PAUSE 

330  PEM  MAKE  MEASUREMENTS  AT  ALL  FREQS  ON  APC-7  SHORT 

340  FOR  K  »  1  TO  N 

330  LET  F3-FCK1-6*F9 

330  FOR  J»1  TO  13 

3  0  0  F  R  E  0  2  t  F  3  > 

320  MEAC 1  ( 1 50 , X , V  I 
340  CPAK i M , V , 2 C J ,KI ) 

330  LET  F3«F3+F3 
3  3  0  N  E  T  J 
7  0 0  NEXT  K 

720  DIRLAV  "CONNECT  APC-7  SLIDING  LOAD" 

740  PAUSE 

750  PEM  MAKE  MEASUREMENTS  AT  ALL  FPEQS  ON  RPC-7  SLIDING  LOAD 
730  FOR  K . 1  TO  N 
7. 3  0  LET  F3.FCI  J-8*F9 
3  0  0  FOP  J « 1  TO  13 
3 2  0  F  R  E  0  2 l F  3 ) 

340  MEAC 1 1 130.X. Y) 

330  CPAK  CI.V.LU.Kl) 

2  3  0  LET  F  2  •  F  3  *  F '? 

3  0  0  r  *  E  X  T  .t 
420  NEXT  I 

>40  DIPLA,  “CONNECT  APC-7  OPEN" 

3  3  0  P  A  U . E 

470  PEN  UAIE  ME  A  2  UPEUENT  I  AT  ALL  FRET'S  ON  APC-7  QREN 


Fig.  13.  Page  1  of  listing  for  Program  1 
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PAGE  2 


■?  S  0 

1  0  0  0 
0 
V 

0 
0 
0 
0 
0 

1  1  .10 
114  0 
1100 
1  1  ?  0 
1  2  0  0 
12  2  0 
'?  ?  >.  0 
0  0  9  0 


r  0  P  K » 1  TO  N 
LET  F3»FCK3-6*F9 
FOP  J.l  TO  13 
P  E  Q  2  (  F  3  ) 

ME  ft  2 1  i  1 50 ,  X , Y ) 

c  pat  ( .  v  ,ot  j  ,k  3 ) 

LET  F3.F3+F9 
next  j 
next  K 

PEM  3  TOR  E  ALL  NECESSARY  LOCAL  VARIABLES 

LET  Ntll-Fl 

LET  NC21-F2 

LET  NI3I.F9 

LET  NC01.S 

LET  NC71.N 

PEM  BRING  IN  NEXT  PROGRAM  AUTOMATICALLY 
CHA  1 1!  i  "  1  I  ANU73  .  S  "  1 


IN  COMMON 


Fig.  14.  Page  2  of  listing  for  Program  1 
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(the  first  measurement  frequency) ,  and  initializes  the  RF 
measurement  equipment.  The  FI's  are  then  calculated  and 
stored.  Lines  500  to  700  perform  measurements  on  the  APC-7 
short  at  all  satelite  frequencies  (frequency  measurements 
around  each  FI  necessary  to  compute  discrete  Double 
Running  Average  for  the  FI's).  Lines  720  to  920  and 
940-1120  perform  the  same  measurements  on  the  APC-7 
sliding  load  and  the  APC-7  open,  respectively.  The  local 
variables  are  then  stored  in  COMMON  to  be  passed  to  subse¬ 
quent  programs,  and  the  program  ends  by  chaining  program 
2  into  memory. 

The  flowgraph  of  program  2  appears  in  Fig.  15  and 
the  listing  in  Fig.  16.  This  program  is  tasked  with 
characterizing  the  open  circuit  and  computing  the  error 
correction  coefficients.  It  begins  by  recalling  the  local 
variables  from  COMMON  (assigning  more  descriptive  names) 
and  by  defining  the  complex  constants  1  +  jO  and  -1  +  j0. 
Then  for  each  frequency,  open  circuit- phase  adjustments 
are  calculated  and  applied  to  the  ideal  reflection  coef¬ 
ficient  of  the  open  circuit  to  obtain  T  in  Eq.  23.  The 
error  correction  coefficients  are  then  calculated. 

Lines  340-400  are  the  code  to  compute  Q  in  Eq.  27,  lines 
420-460  to  compute  e^  in  Eq.  29,  and  lines  480-520  to 

in  Eq.  30.  These  coefficients  then  replace  the 


compute  e 


1  REM  THU  PROGRAM,  DEVELOPED  ON  U/^.89,  ID  DESIGNED 

3  REM  TO  CALCULATE  THE  ERPOR  CORRECTION  COEFFICIENTS 
5  REM  FOP  THE  A  P  C  -  7  CALIBRATION 

4  REM  REV  1  1  '  7  <-  8  O  REM  3  '  1  3 ' 8  1  1  :  A  N  U  7  3  .  S 

3  0  COM  *[13,?l),LU3,51)t0CI3,!ll 

40  COM  F  C  5  1  3  ,  N  [  1  O  3 

9  0  REM  RECALL  LOCAL  VARIABLES  FROM  COMMON 
100  LET  F9.N C 3  3 
130  LET  3 . N I  6 3 
14  0  LET  N  .  N  C  7  3 

150  REM  DEFINE  COMPLEX  CONSTANTS  (1,0)  AND  (-1,0) 

1 00  CPAK (1  , 0  ,  D2 ) 

180  CPAk(-1.0.D3) 

190  REM  GENERATE  ERROR  CORRECTION  COEFFICIENTS  FOR  ALL  FREQ 

3  0  0  FOR  K  «  1  TO  N 

330  LET  F3-FCK3-6*F9 

340  FOR  J»t  TO  13 

350  REM  CORRECT  APC-7  OPEN  CIRCUIT  PHASE 

300  LET  C9-5.02E-05+1.126E~l4*F3t2 

330  LET  C9.-C9 

300  PCFT  i 0 , F  3 ,C9  ,  D2  ,P1 ) 

310  REM  COMPUTE  ERROR  COEFFICIENTS 

311  REM  EO.EO0  ,  E 1 -E  1  1  ,E2-E2 1 
320  LET  E0«LCJ,K3 

340  CSUB ( SC  J  ,K 3  ,E0  ,N1 ) 

360  CSUB ( OC  J  ,K  3  ,  E  0  ,  D I  ) 

3  30  CDIVlNl  ,D1  ,N1) 

400  CMPV (PI ,N1 ,Q5) 

430  CADD (Q5  ,D2  ,N1  ) 

440  CSUB(P1, 05,01) 

460  CDIVlNl, D1.E2) 

4  30  CSUB ( E0  ,  S C J  ,  K 3  ,N 1 ) 

500  CADD < D2  ,  E2  ,  D 1  ) 

520  CMPV (Ml  .  D 1  ,  E 1  ) 

531  PEM  SAVE  ERROR  CORRECTION  COEFICIENTS  IN  COMMON 

540  LET  SIJ,K]»E0 

560  LET  Lt  J  ,K  3-El 

580  LET  u  C  J  ,  K  I  .  E  2 

590  LET  F3.F3«-F9 

600  NEXT  J 
620  NEXT  K 

9930  REM  LOAD  NEXT  PROGRAM  AUTOMATICALLY 
9  9  9  0  C  H  A  I  N (  "  1  :  A  N  U  ~  4  .  :  "  i 


Fig.  16.  Listing  for  Program  2 
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APC-7  measurements  in  COMMON  since  they  are  no  longer 
needed,  and  program  3  is  chained  into  memory. 

Figures  17  and  18-19  display  the  flowgraph  and 
listing  for  program  3  which  makes  system-corrected 
measurements  on  adapter,  reference  line,  and  load  and 
computes  the  Double  Running  Average  for  each  FI.  It 
starts  with  virtually  the  same  variable  recalling  and 
initialization  as  program  2.  The  adapter,  load,  and  line 
are  then  measured  and  corrected  for  system  errors  (lines 
520-580  represent  Eq.  31) .  Lines  680-960  represent  the 
code  for  computing  the  first  running  average  sum.  Note 
the  application  of  the  sliding  running  average  in  lines 
900-960.  The  second  running  average  sum  is  calculated  in 
lines  980-1100  with  the  average  being  completed  in  lines 
1140-1160.  COMMON  is  then  rearranged  so  that  a  transition 
to  single  subscripted  arrays  can  be  accomplished  and  so 
that  a  large  block  of  memory  can  be  freed  in  program  4 
since  only  information  at  the  FI's  is  required  from  this 
point  on.  Program  4  is  now  chained  into  memory. 

The  flowgraph  and  listing  of  program  4  are  displayed 
in  Figs.  20  and  21,  respectively.  This  program  initi¬ 
alizes  the  equipment,  measures  the  adapter  and  short, 
corrects  the  data  for  system  error,  and  chains  program  5 
into  memory.  At  this  point  all  the  information  necessary 
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Fig.  17 


Flowchart  for  Program  3 
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PAGE  1 

1  PEM  T  H  1 1  P P 0 G ROM  ,  D E  V ELOPED  ON  11/5/80,  IS  DESIGNED 

2  PEN  TO  COMPUTE  THE  DOUBLE  PUNNING  AVERAGE  FOP  THE 
3-  P  E  M  ADOPTER  CORRECT  I  0  M 

4  P E  M  REV  11'?'80  PEM  3/13/81  1  :  A  N  U  7  4  .  S 


2  0 

COM  A  C 1 3 ]  , 

BCI3J,CC133,DC13 

3  ,  E  C  1  3  3 

4  0 

COM  2  C  1  3 , 5 

1]  ,  L  C  1  3 , 5  1  3  ,  0  C  1  3 

,513 

6  0 

C OM  F  C  5 1  ]  , 

N  C  1  0  ] 

3  o 

P E  M  RECALL 

LOCAL  VARIABLES 

FROM  COMMON 

1  0  0 

LET  F 1  *  N  C 

1  3 

1  2  0 

LET  F  2  «  N  C 

2  3 

1  4  0 

LET  F  9 ■ N  C 

3  3 

1  6  0 

LET  N«N  C 7 

3 

17  0  PEN  INITIALIZE  THE  MEASUREMENT  E  0 IJ  I  P  M  E  H  T 
130  FCALF<FC13-**r9« 

200  BCNT1  ( F  C  1  ]-o*F9 J 
2  2  0  : S  E  L 1  '  1  n 

240  UA I T  (50) 

250  PEM  CLEAR  DISPLAY 

2  6  0  CLEAR  1  O' 

270  REM  DEFINE  COMPLEX  CONSTANT  (1,0) 

230  C  P  A K  t  1  ,0 , 02 ) 

260  DSPLAY  "CONNECT  SMA  LOAD  AND  LINE" 

380  PAUSE 

400  FOR  K  *  1  TO  N 

420  LET  F3-FCK]-6*F9 

430  REM  MAKE  TWO  CYCLE  MEASUREMENTS  FOR  OUTPUT  POINT  K 

431  REM  ON  ADAPTER,  LINE,  AND  LOAD 
4  4  0  FOR  J  »  1  TO  13 

46  0  F R E Q 2 i F  3 ) 

4  3  0  M E  A 1 1  l  1 5  0  ,  X  ,  Y ) 

5  0  0  C  P  A  K  (  X  ,  Y  ,  Z  ) 

510  PEM  PERFORM  ERROR  CORRECTION  ON  ADAPTER,  LOAD,  AND  LINE 
520  C3UBiZ,SCJ,K3,Dl) 

5  4  0  C  D  I  V  i  L  C  J  ,  K ]  ,  D 1  ,  D  1) 

560  C ADD'D!  ,  0  C  J  ,  K  J , D 1  ) 

530  CDIViD2,Bl,ACJ]) 

600  LET  BCJ3«REA(ACJ3) 

6  20  LET  C  C  J  3  -  I  M  G  (  A  C  J  3  ) 

640  LET  F3-F3+F9 

660  NEXT  J 

670  REM  I ET  UP  DOUBLE  RUNNING  AVERAGE  FOR  POINT  K 
6  7  1  PEM  MAKE  SEVEN  M E  0  S  U  R  E  M E  N T  S  PER  CYCLE 
6 s o  let  M * 3 

6  9  0  PEM  COMPUTE  FIRST  PUN  A  V  G  SUM  FOR  POINTS  NEEDED  TO  COMPUTE 
■=■  l  PEM  SECOND  PUNNING  AVERAGE  FOR  OUTPUT  POINT  K 

7  0  0  l  E  T  1=M  +  1 

7  2  0  LET  DCM«£CI-M]'2 
7  4  0  LET  E  C  I  3  «  C  C  I  -  M  3  *  2 
7  6  0  FOR  .*  %  I  -  M  «•  1  TO  I  ♦  M  -  1 

7  0  LET  D  [  I  J  =  D  t  J  3  *  B  [  J  J 

3  0  0  LET  EtIJ=ECIl«*CtJ] 

3  2  0  LET  J  1  * J 
3  4  0  NEXT  J 

3  6  0  LET  DCIJsDCIl+PtJlMl/S 
3:0  L£~  £  :  I  1  * E  C  I  1  «■  C  i  J  1  ♦  1  3  •  2 

-  -  m  '  Z.\  CONFUTE  SLIDING  FUN  A  G  TUtIC  FOP  PEST  OF  POINTS  DEEDED 

j  1  F  [  "  ■  *1  c  M  F  U  T  E  :  c  *  o  N  D  -  IJN  M  V  F  M  P  P  0  I  N  T  K 

r.-.i  2  TO  I  '  -  M 


Fig.  18 


Page  1  of  listing  for  Program  3 
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F  Cm;  E  2 

■=•  J.  LET  D  [  I  3  *  D  C  I  -  1  3  -  B  C  I  -  M  -  1  3  -  E  C  I  -  N  3  B  C  I  -  M  -  1  J  1  J  *  B  [  I  -  M  ■  1  ;  ■  2 

■?  -i  let  ec  i  3  *ec  i  -  i  i  *  i  -  c  c  i  -  i  j-cc  i-«  J*c  c  i-«-i  ♦  jj  j  *c  c  i  -  ■  ; 

?  4  0  NEXT  I 

1  PE  M  C  0  II P  U  T  E  EE  C  0  N  D  RUN  A  V  G  CUM  EuR  POINT  K 

9  <3  0  LET  I  =  £  »  ri  +  1 

10  0  0  LET  U  =  D  C  I  -  M  3 ' 2 

10  2  0  LET  T  .  E  t  I  -  M  3  •'  2 

1  0  4  0  POP  J  =  I  -  M  ♦  1  TO  IP  IT  -  1 

1  0  p.  0  LET  U  *  U  *  D  C  -I  I 

1000  LET  T«T*ECJ3 

1  1 0  0  LET  j l , J 

1120  ne::t  j 

112  1  PEP  D  I  IDE  BY  TOTAL  ME  I GH  TED  BUM  TO  GET  SECOND  PUN 
112  2  P  E  IT  A  V  G  POP  POINT  K 

1.40  lET  U»iU  +  DEJ1*13'2i/’36 

11-0  LET  T*iT+ECJ1+13- 21/36 

1  1  4.  1  PEN  :  T  o  P  E  PUN  A  V  G  LOAD  AND  R  E  A  P  P  ft  N  G  E  COMMON  IN  PREPARATION 

11-2  PEN  POP  VARIABLE  REASSIGNMENT:  IN  NEXT  PROGRAM 

1  i  %  0  C  P  A  K  I  IJ  ,  T  .  0  t  1  2  ,  K  ]  I 

1  2  0  0  -  E  T  0  C  3  ,  l:  3  *  S  C  7  .  K  3 

1  2  2  0  -El  0  C  9  .  I  3  =  L  C  7  .  M 

12  40  LET  0 C 1 0 . 1  3  =0  C  r  .  I  3 

12G0  .-EOT  k 

9  9  SO  PE  IT  LOAD  NEXT  PROGRAM  AUTOMATICALLY 
'?  9  9  0  0  H  A  I  N  l  "  1  :  A  N  U  7  5  .  S  "  i 


Fiq.  n 


Pnon  2  of  listinq  for  Program  3 


PEM  THIS  PROGRAM  ,  DEVELOPED  ON  11/^  30,  IS  DESIGNED 
P E  M  TO  COMPUTE  THE  ADAPTER  CORRECTION  COEFFICIENTS 
REM  REV  1  1  '  2  4  -  3  O  P  E  M  3 ' 1  3 ' 3  1  1  :  A  N  U  7  3 . 3 

2  0  COM  C  C  2 ♦ 5 1  ]  ,  S C  5  1  3  ,  L  C  5  1  3  ,  MC  5  1  3 
4  0  COM  F  C  5  I  3  , N  C  *  0  3 

41  REM  RECALL  LOCAL  VARIABLES  FROM  COMMON 
100  L  E  T  N*N [ 7 3 

MO  REM  INITIALISE  MEASUREMENT  EQUIPMENT 
12  0  F  C A  L  F  i  F Z  1  3  > 

140  ECHTUrtlD 
1*0  3 E  L  1  I  1  1  -I 

1  3  0  M  A  I  T  t  5  0  :» 

2  0  0  D : P L  A  Y  "CONNECT  SNA  SHORT" 

2  2  0  P  A  U  S  E 

2  3  0  REM  DEFINE  COMPLEX  CONSTANT  (  1  , 0 1 
2  4  0  CPAK  i  1  . 0  .  D 2 ) 

250  PEM  MAKE  CORRECTED  MEASUREMENTS  ON  ADAPTER  AND  SHORT 
2 6 O  FOR  I  =:  l  TO  N 

2  3©  F  R  E  Q  2  l  F  C  I  3  > 

3  00  ME AS  1  i  1 5  0  ,  X  ,  Y ) 

32  0  CPAK  (  X  ,  V  .  2  :■ 

3  4  0  C S  U  B • .  Z  ,  C  C  1  ,  I]  ,  D  1  ) 

3*0  CDIVICC2.I3.D1.D1* 

3  8  0  CADDID1  ,  C  C  3  ,  I  3  ,D1  ) 

4O0  C  D I  V ( D  2  ,  D  1  .  S  C  I  3  i 

4  4  0  n  e : :  T  i 

•?  ?  3  0  PEM  LOAD  N  E  : ;  T  PROGRAM  AUTOMATICALLY 
•?  ?  ’?  0  C  H  A  I  N  1  "  1  ‘A  H  Ll  7  *  .  S  "  > 


Fig.  21.  Listing  for  Program  4 
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to  make  fully  corrected  measurements  has  been  obtained, 
and  the  system  is  prepared  to  measure  the  complex  re¬ 
flection  coefficient  of  DUT's. 

Figures  22  and  23  show  the  flowgraph  and  listing 

for  program  5  which  has  the  job  of  measuring  a  DUT  and 

fully  correcting  for  system  error  and  the  adapter.  After 

initializing  it  makes  measurements  on  a  DUT  and  corrects 

it  for  system  error  in  the  usual  way.  Lines  460-540 

computes  the  value  of  El  in  Eq.  9.  Lines  553-556  normal- 

* ■  i  2  0 

ize  El  to  make  sure  that  it  is  indeed  equal  to  e  J 

Lines  560-640  complete  the  calculations  of  Eq.  10  yielding 

the  fully  corrected  reflection  coefficient  of  the  DUT. 

This  result  will  now  be  output  in  program  6. 

Program  6,  whose  flowchart  and  listing  appears  in 
Figs.  24  and  25-26,  provide  plots  of  Real  (T)  vs.  freq. , 
Imag  (D  vs.  frequency,  Mag  (T)  vs.  frequency,  and  Ang  (D 
vs.  frequency.  It  will  then  provide,  if  desired,  a  com¬ 
plete  printout  of  the  above  mentioned  parameters. 

Program  5  is  then  chained  into  memory  and  the  system  is 
ready  to  measure  another  DUT  of  the  same  connector  type. 
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PAG  E  1 

1  t'DI  THIS  PROGRAM,  DEVELOPED  OH  11  5 --SO.  IS  DESIGNED 

2  REM  TO  ADAPTER  CORRECT  THE  MEASUREMENT  OF  THE  UNKNOWN 

3  REM  REV  1  2  ■  2  2  '  3  0  PEN  3M3/3I  1  i  A  N  u  7  6  .  S 

20  COM  CC3,513,S£513,LCT13.MC513 

40  COM  F  C  5  1  3  .  N  [  l  0  3 

50  REM  RECALL  LOCAL  VARIABLES  FROM  COMMON 
1O0  LET  N  •  N  C  7  3 

110  REM  INITIALIZE  MEASUREMENT  EQUIPMENT 
ISO  F  CALF  (  F  C  1  3  I 
14  0  B  C  N  T  1  (  F  C  1  3  I 
160  SSEL  1  Ul) 

ISO  WAIT  i50i 

185  REM  CLEAR  DISPLAY 

190  CLEAR ( O I 

200  D'PLAV  "CONNECT  SMA  UNKNOWN* 

220  PAUSE 

230  PEM  DEFINE  COMPLEX  CONSTANT  (1,0) 

24C  CP  AK  !  1  ,  O  ,  D2  ) 

250  REM  MAKE  ERROR  CORRECTED  MEASUREMENTS  ON  UNKNOWN 
260  FOR  1-1  TO  N 
280  FREQS ( F I  I  3  I 
300  MEAS 1  ( 1 50  .X  ,  Y) 

320  CPAK (X ,Y ,2) 

340  CSUB iZ,CCl,I3,Dl) 

360  CDIV (CC2  ,  I  3  ,D1  ,D1 ) 

380  CADD ( D 1  , C  C  3  ,  I  3  ,  D 1  I 
400  CD  I  V  (  D2  ,tl  1  ,MC  I  3  ) 

420  NEXT  I 

43C  PEM  MAKE  ADAPTER  CORRECTIONS  FOR  ALL  OUTPUT  POINTS 

440  FOR  I»1  TO  N 

460  CSUB (LC I  3  ,St  I  3  ,D1 ) 

430  CPAKCREAILCI3)  ,-IMG (LC  I  3  3  ,N2> 

500  C M P Y (SC  I  3  ,N2  ,N1 ) 

520  CSUB ( D  2  ,  N 1  ,N1 I 
540  CDIViNl , D 1 ,N3) 

553  CPAK  (MAG (N3)  ,0  ,N4) 

556  C  D I  V l N3 , N4 , N3 I 
5  6  0  C  S  U B  < M  C I  3  ,  L  C  I  3  ,  N  1  I 

5  3  0  C  M  P  Y  !  M  C  I  3  , N  2  ,  D 1  I 
6 ■'  0  CSUB  i  D2  ,  D  1  ,  D  1  ) 

620  CD  I  V (N1  , D 1  ,N  1  I 

6  4  0  C  M  P  Y i N  1  , N  3 , M  C  I  3  > 

660  NEXT  I 

9930  PEM  LOAD  NEXT  PROGRAM  AUTOMATICALLY 
9990  CHAIN'  "1  ! A  N  U  7  7 . S " 1 


Fig.  23.  Listing  for  Program  5 
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PAGE  l 

1  PCM  THIS  PROGRAM.  DEVELOPED  ON  ll/-5'80,  IS  DESIGNED  ! 

3  REM  TO  OUTPUT  THE  CORRECTED  MEASUREMENTS 

3  REM  REV  11.'?.  30  REM  12x22x30  liANU?7.S  1 

20  COM  CC3, 513, SC513,L£513,MC 313, FC513, Nil  03 

30  REM  RECALL  LOCAL  VARIABLES  FROM  COMMON 

100  LET  F 1 >N C  1  3 

120  LET  F2«N  t  2  3 

140  LET  $  >  N  [  6  3 

100  LET  N.NC73 

170  REM  INITIALIZE  GRAPHICS  DISPLAY  AND  BUFFER 
130  BUF l 1 00 ) 

2  08  clear  i  0  j 

2  20  DCPLAY  "LI  HAT  ARE  BOTTOM  AND  TOP  LIMITS"; 

240  BELL 

2  3  0  INPUT  B 1  , T  1 

230  LET  S1«(T1-B1)/-10 

230  REM  CLEAR  DISPLAY,  SCALE,  AND  PLOT  REAL  PART  OF  REFLECTION 
300  C  L  E  A  R ( 0  j 

3  2  0  S  C  A  L  E ( F  1  , F  2  .  B  1  ,  T  n 
340  LET  F3. ( F2-F 1 ) x  10 
300  S AXES  1 F3 , 3  1  ) 

330  DSPLAY  VTAB (  1 5 )  ,  T AB ( 0 )  , “ REFL" 

4O0  DSPLAY  VTAB  I  3 1 )  .  TAB ( 23 )  . “ FREQ (MHZ )  » 

420  LABEL  VTAB ( 16)  , TAB (0)  , “REAL" 

440  BLOCK  l  B2 ) 

460  FOR  1«1  TO  N 

430  PLOT (Ft  I  3  ,REA (Mt I  3 1  ,2) 

5O0  NEXT  I  ! 

520  PAUSE 

530  REM  CLEAR  REAL  PLOT  AND  PLOT  IMAGINARY 
540  CLEAR i 2 ) 

500  SAXES ( F3  .  S  1  )  i 

530  LABEL  V T A B (  1  6 )  ,  T A B ( 0 )  ,  " I M AG " 

600  BL0CK(B2) 

620  FOR  I«1  TO  N 

640  PLO f  i F C  I  3  , I  MG ( M C I  3  )  ,  2 ) 

660  NEXT  I 
630  PAUSE 

700  CLEAR (01  ! 

720  DSPlAY  "UHAT  ARE  BOTTOM  AND  TOP  LIMITS  FOR  MAG";  ! 

740  BELL 
760  INPUT  B2.T2 

770  REM  CLEAR  DISPLAY  AND  PLOT  REFLECTION  MAGNITUDE 

7  3  0  C  L  E  A  R  (  0  i  ; 

';00  LET  ;2«iT2-B2)/lG  ; 

320ZCftLE«Fl,F2,B2,T2)  ! 

3  4  O  Z  ft  ^  E  2  (  F  3 , 1 2  >  t 

3*0  Ii  *  P  t  ft  7  7  T  Q  B  i  1  5  )  ,  T  ft  B  (  0  )  ,  "  R  E  F  L  " 

3  30  D  Z  P  L  ft  Y  V  T  ft  B  i  3  1  )  ,TftB(28)  ,  M  FREQtMHZ)" 

9  O  O  LftBEL  V  T  ft  B  <  l  6  J  ,  T  ft  B  f  Q  )  ,  " MAC" 

92  0  BLOCK  f  B 3 ) 

940  POP  I.l  TO  M 
9  *  O  P  L  0  T  «  F  :  I  ]  ,  M  ft  0  i  M  C  I  3)  ,  2  > 

9  3  0  N  E , :  T  I 

1  0  0  0  "  ft  U  I  E 

10  1  9  c  E  ’*  iIlEuP  rift  0  5!  I  T'jrE  PLOT  PNC  PLOT  REFLECTION  ft  NC  L  E  ( 

1  o  2  0  C  L  E  ft  P  '•  2  i  I 

1  0  4  0  l[i  t  '  =  -  l  •?  0 


Fig.  25.  Page  1  of  listing  for  Program  6 
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P  P  5  E  2 


1 k"'  6  0 

LETT  T3-1S6 

1  0  -3  0 

LET  53»<T3-B3)--10 

1  1  OO 

SCALE  CF1  ,F2 .E3  ,T3> 

112  0 

JAKES i F3 . S3) 

1140 

LABEL  VTABtl6),TAB(0t,"ANG" 

l  160 

BLOCK  i  B2 ) 

l 

1130 

FOR  I » 1  TO  N 

1  2O0 

PLOT (  F  C  I  3  , ANG (MC I  3  )  ,2) 

1 

1  2  2  0 

NE::T  I 

1 

12  30 

PEM  CLEAR  DISPLAY  AND  PRINT  ALL  RESULTS  FOR  ALL 

123  1 

REM  FREQS  IF  REQUESTED 

1 

124  0 

PAUSE 

20  0  0 

C  L  E  A  P  (  O  ) 

J 

2  0  20 

DC  PLAY  "IJANT  TO  PRINT  RESULTS(1-YES,2-N0)”j 

2040 

BELL 

20  6  0 

INPUT  HI 

2080 

IF  Hl»l  GOTO  9000 

209O 

REM  INITIALISE  AUTOMATIC  GRAPHICS  BUFFER  PACE  TURNER 

2  100 

TRAP  8  GOiUB  3000 

2120 

CLEAR  ( 0  ( 

2140 

BUF  (3 00  J 

REFL(WBG) " |  I 

2180 

DSPLAY  "FREQ 'MHZ)  REFL(REA)  REFL(IMC) 

220O 

DSPLPY  "  REFL(ANG)” 

2220 

DSPLAY 

1 

2240 

FOR  I»1  TO  N 

r 

2250 

REM  SET  OUTPUT  FORMAT  IN  BASIC  NOTATION 

1 

226  0 

FDSPtFI I 3 , ? ,0) 

1 

2280 

DSPLAY  '•  "  ; 

2300 

FDSP I  REA ( MC I  3 )  ,6,3) 

2320 

DSPLAY  "  "! 

2340 

FDSP (  IMG (MC I 3  ) ,6,3) 

2360 

DSPLAY  "  "  } 

2380 

FDSP ( MAG ( M  C I 3 )  ,6,3) 

2400 

DSPLAY  '•  “I 

! 

2420 

FDSP  l  ANG (MCI))  ,6,1) 

I 

2440 

DSPLAY 

1 

2460 

NEXT  I 

1 

2470 

COPY 

1 

24  75 

PAGE 

2  4  7  7 

PAGE 

1 

2480 

GOTO  9000 

2  99  0 

REM  SUBROUTINE  FOR  AUTOMATIC  PAGE  TURNINC  WHEN 

GRAPHICS  ! 

2  9  9  1 

PEM  BUFFER  BECOMES  FULL 

3  0  0  0 

C  0  P  Y 

i 

3  0  1  0 

PAGE 

2  015 

PAGE 

2  0  2  0 

C  L  E  A  P  (  0  ) 

2  0  4  0 

fetij  r  h 

2  '■*  0 

PEM  load  measurement  ppogpam  automatically 

j 

9  0  0  0 

C  H  A  1 N  '  “  I  :  A  N  U  T  15  .  S  *  ' 

Fig.  26.  Page  2  of  listing  for  Program  6. 
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CALIBRATION  PROCEDURES 

The  calibration  procedures  are  not  difficult  to 
accomplish.  In  fact,  the  computer  gives  rather  explicit 
instructions  for  each  step,  and  always  signals  with  a  bell 
when  operator  action  is  required. 

Only  the  first  program  must  be  loaded  into  memory 
manually  (the  rest  are  sequenced  automatically) .  It  first 
asks  for  the  electrical  length  of  the  reference  line  being 
used  for  the  load  characterization.  It  then  returns  the 
minimum  and  maximum  frequency  limits  and  asks  for  the 
start,  stop,  and  step  frequencies.  This  information  is 
then  used  to  calculate  the  FI’s.  The  operator  is  next 
prompted  to  connect  an  APC-7  short,  an  APC-7  sliding 
load,  and  then  an  APC-7  open.  After  calculating  the 
system  error  correction  coefficients,  the  program  asks 
that  the  adapter,  reference  line,  and  load  be  attached 
after  which  it  computes  the  Double  Running  Averaged  load 
characterization.  Connection  of  the  adapter  and  short  is 
then  requested,  and  the  calibration  is  complete  and  the 
system  is  prepared  to  make  fully  corrected  measurements 
on  any  DUT. 

After  attaching  a  DUT  to  the  adapter,  the  system 
makes  measurements,  corrects  these  measurements,  and  then 
asks  for  the  bottom  and  top  limits  of  the  r(Real)  and 
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r(Imag)  plots  versus  frequency.  After  displaying  these 
plots,  the  program  requests  the  bottom  and  top  limits  for 
the  r (Mag)  plot  versus  frequency.  This  is  output  and  the 
F (Ang)  plot  versus  frequency  automatically  follows. 

The  operator  can  now  choose  if  a  printed  listing 
of  the  frequencies  and  above  parameters  are  to  be  output. 
Whether  or  not  the  listings  are  requested,  the  computer 
automatically  asks  for  another  DUT  to  be  attached  to  the 
adapter,  and  the  sequence  repeats  until  the  operator 
interrupts  the  program  by  pressing  the  BREAK  button  on  the 
console. 
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COMPARISON  OF  TRANSMISSION  LINES 

As  was  stated  earlier,  use  of  the  Double  Running 
Average  for  characterization  of  an  ideal  load,  for  the 
first  time,  permits  direct  comparison  of  transmission 
lines  and  connection  systems.  For  instance,  two  or  more 
lengths  of  transmission  line  can  be  directly  compared 
against  each  other  to  obtain  an  exact  impedance  match  at 
the  frequency  or  band  of  interest. 

However,  this  procedure  cannot  be  conducted  using 
the  calibration  and  measurement  programs  as  already  de¬ 
scribed.  They  must  be  modified  and  these  modif ications 
will  now  be  presented.  Programs  1  and  2  remain  completely 
unchanged  and  program  3  requires  only  a  change  in  the 
CHAIN  statement  to  the  new  program  4  (Program  4A) . 

Figures  27  and  28  show  the  flowgraph  and  listing 
for  program  4A.  This  program  recalls  the  local  variables 
from  COMMON,  opens  a  disc  file  for  data  storage,  and  then 
asks  for  identification  of  the  transmission  line  for  which 
data  is  being  stored  (up  to  4  lines  can  be  compared  at  one 
time  with  this  program) .  The  line  data  is  then  stored  on 
the  disc  along  with  the  error  correction  coef f icients . 

The  program  then  halts.  If  another  line  is  to  be  measured 


Fig.  27 


Flowchart  for  Program  4A 
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4  0 


r EM  T Hi:  FROG PPM.  DEVELOPED  ON  il/r/80,  I?  DET 
P  E  M  J  T  0  P  L  THE  PVERPC-ED  LINE  DP  T  D  OH  D  I  Z  C  F  0  P  C 
p E  M  OF  DIFFERENT  LINE  3  D I R  E  C  T L V 

R  EM  PE  V  i  1  2  4*-  3  1  P  E  M  3  '  X  3  3  1  1  :  Q  N  U  7  3 

PE  M  THi:  P  P  0 6  R  P  M  MUST  REPLACE  1  :flNU75  .C  W  HEN 
REM  GENERPTING  THE  DP  TP  FILE  FOR  COMPflR I  CON  OF 
P E  H  LOmDS  PHD  LINES 
•:  0  M  C  :  3  «  5  1  3  ,  3  C  5  1  ]  ,LC513 
•:  0  M  F  I  5  1  3  ,  N  I  l  0  ] 

P  E  M  ? E  C P  L  L  L  0  C  ft  L  VPRIPBLES  FROM  COMMON 


■SHED  T 

M  P  P  P  I  : 


100  LE.T  N  *  h  C  7  3 

I  JO  FEM  OPEN  DISC  FILE  FOP  DPTQ  STORPGE 

:  ^  0  0  PEN*  *  0  ,  "  1  :  P  N  U  7  *  .  D  "  ,  E  1 

1  : 0  p  E  M  D  P  T  P  F  I j P  or  ~0  FOUR  LINES  C  P  N  F  E  STORED 

.  4  0  D : R  L  P  Y  "WHICH  LINE  *  1 »  0  U  P  L  ,2- SHORT, 3-LONG  ,4*. SPLINE) 

I  *s  0  FELL 

1  :?  0  I  N  P  0  T  0 

1  •?  o  p  e  m  d  e  p  o  :  I  t  d  p  t  p  in  Dire  rue 

2  0  0  D  P ’  ' E  i  1 0  ,  0-1  j  *  3 1 ♦ 1  .Li  1  1  ,  j  l  ,E' 

2  2  0  D  P  I  T  E  •  10,2  0  ,  C  C  1  .1  3  ,  1  5  3  ,  E  ' 

5  =»  ■?  0  END 


N 


j  Fin.  28.  Listinq  for  Proqr  4A. 
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then  program  1  must  be  reloaded  into  memory  and  the 
sequence  repeated.  After  all  lines  have  been  measured, 
program  SA  can  be  loaded  and  run. 

The  flowgraph  and  listing  of  this  program  are  dis¬ 
play  in  Figs.  29  and  30.  This  program  must  measure  the 
adapter  and  short,  correct  it  for  system  errors  and  store 
this  data,  the  frequency  list  and  COMMON  variables  on  the 
disc.  All  the  data  necessary  for  comparing  the  lines  has 
now  been  stored  on  disc. 

Programs  4B,  5B  and  6  are  used  to  perform  the 
actual  line  comparisons  and  to  output  the  results.  The 
sequence  is  started  by  loading  and  running  program  4B. 

The  program  sequence  progresses  automatically  from  this 
stage.  Program  4B  [Figs.  31  and  32]  begins  by  loading  all 
the  data  (except  line  data)  into  COMMON.  It  then  asks 
which  line  will  be  used  as  the  reference  line  and  loads 
the  appropriate  data  into  COMMON. 

Program  5B  (Figs.  33  and  34)  is  then  automatically 
chained  into  memory.  This  program  requests  which  line  will 
be  used  as  the  DUT  and  loads  that  data  into  COMMON  from 
the  disc.  It  then  performs  the  adapter  correction,  which, 
in  this  case,  actually  represents  the  line  comparison. 
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PAGE  1 

REM  THIS  PROGRAM  ,  DEVELOPED  ON  ll/'S/'Se,  IS  DESIGNED 
REM  TO  MEASURE  AND  STORE  DATA  FOR  LINE  COMPARISON 
REM  REV  U''24/80  REM  S^'lS-'Sl  1|ANU?9.3 

REM  THIS  PROGRAM  MUST  REPLACE  1iANU75.S  WHEN 
REM  STORING  THE  ADAPTED  SHORT  AND  FREQ  DATA  FOR  AVERAGED 
REM  LOAD  AND  LINE  COMPARISON 
COM  CCS, 51  I  ,SC51 3  ,LC51 3  ,MC51I 
COM  FC51 3 ,NC 103 

REM  RECALL  LOCAL  VARIABLES  FROM  COMMON 
3  LET  N.NC7I 

3  REM  INITIALIZE  RF  EQUIPMENT 
3  FCALF(FCII) 

3  BCNT1  (Ft  1 3 ) 

3  S3EL1 (11) 

3  WAIT  (50) 

3  REM  MEASURE  AND  STORE  ADAPTER  AND  SHORT 
3  DSPLAY  "CONNECT  SMA  SHORT" 

3  PAUSE 

3  REM  DEFINE  COMPLEX  CONSTANT  (1,0) 

3  CP AK (1,0, D2 ) 

3  FOR  I«1  TO  N 
3  FREQ2(FCI3) 

3  ME AS  1  (1 50 , X , Y ! 

3  CPAK (X,Y ,2) 

3  REM  CORRECT  MEASUREMENTS  FOR  SYSTEM  ERRORS 
3  CSUB (2 ,CC 1 , I S ,01) 

3  CDIV(CC2,I3  , D 1  , D  1  ) 

3  CADD (D1  ,  C  C  3  ,  I  1  ,  D  1  ) 

3  C  D  I  V  (  D2  ,  D  1  ,  S  C  I  I  ) 

3  NEXT  I 

3  REM  OPEN  DISC  FILE 
3  OPEN ( 10, " 1  »ANU?3  .D" ,E) 

3  REM  STORE  SHORT,  FREQUENCY , AND  COMMON  DATA  ON  DISC 
3  DRITE ( 10 ,353 , SC  1  I  ,51  ,E) 

3  DRITE ( 10 ,409 ,FC 1 1 ,61  ,E) 

90  END 


Fig.  30.  Listing  for  Program  5A. 
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PAGE  1 

1  REM  THIS  PROGRAM,  DEVELOPED  ON  11^5'S0,  IS  DESIGNED 

2  REM  TO  SET  UP  REFERENCE  LINE  FOR  LINE  COMPARISON 

3  REM  REV  5^8^81  REM  3^13^81  1IANU80.S 

4  REM  THIS  PROGRAM  MUST  REPLACE  liANU75.S  WHEN  COMPARING 

5  REM  ONE  AVERAGED  LOAD  AND  LINE  TO  ANOTHER  AFTER  THE 

6  REM  DATA  FILES  HAVE  BEEN  GENERATED 
20  COM  CC3 ,51 1 ,SC51 I ,LC31 3 ,MCS1 1 

40  COM  FC31 3 ,Nt 103 

50  REM  OPEN  DISC  FILE  WHERE  DATA  IS  STORED 
60  OPEN ( 10 , " 1  I ANU73 . D"  ,E ) 

70  REM  RECALL  SHORT  DATA  FROM  DISC 
30  DREAD ( 10 ,338 ,SC  13  ,51  ,E ) 

35  REM  RECALL  FREQUENCY  DATA  AND  COMMON  VARIABLES  FROM  DISC 
30  DREAD (  10 ,409 ,FC 1 J  ,61 ,E) 

82  REM  RECALL  ERROR  CORRECTION  COEFFICIENTS  FROM  DISC 
93  DREADt 10,203,CC1 ,13 ,133,E> 

97  REM  RECALL  LOCAL  VARIABLES  FROM  COMMON 
100  LET  N.NC7J 

140  DSPLAY  “WHICH  LINE  TO  BE  USED  AS  REFERENCE* 

160  DSPLAY  " ( 1-OURS ,2-SHORT ,3-LONG ,4-SPL INE) * | 

180  BELL 
200  INPUT  Q 

210  REM  READ  IN  REFERENCE  LINE  ORTA 
220  DREAD <10,<Q-1)*31+1,LC1I,31,E) 

9000  REM  AUTOMATICALLY  LOAD  NEXT  PROGRAM 
9990  CHAIN ( “ 1  i ANU81  .S* 1 


Fig.  32.  Listing  for  Program  4B 
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Pig.  33.  Flowchart  for  Program  5B 
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PACE  1 

1  REM  THIS  PROCRAM,  DEVELOPED  ON  11/5/80,  IS  DESIGNED 

2  REM  TO  COMPLETE  LINE  COMPARISON 

3  REM  REV  5/8/81  REM  3/13/81  I1ANU8I.S 

4  REM  THIS  PROGRAM  MUST  REPLACE  11ANU76.S  WHEN  COMPARING 

5  REM  ONE  AVERAGED  LOAD  AND  LINE  TO  ANOTHER  AFTER  DATA 

6  REM  FILE  GENERATION 

28  COM  CC3.51 3  ,SC51 I  ,LC51 I  ,MC51 I 
40  COM  FC513.NC103 

98  REM  RECALL  LOCAL  VARIABLES  FROM  COMMON 

100  LET  N-NC73 

180  REM  CLEAR  DISPLAY 

190  CLEAR ( 0 ) 

195  REM  OPEN  DISC  FILE  WHERE  LINE  DATA  IS  STORED 
200  OPEN ( 10 , " 1 iANU7S .D" ,E> 

220  DSPLAY  “WHICH  LINE  TO  BE  USED  AS  UNKNOWN* 

240  DSPLAY  “( 1-OURS ,2-SHORT  ,3-LONG , 4-SPLINE) 

260  BELL 
280  INPUT  0 

290  REM  READ  IN  DATA  FOR  OUT  LINE  FROM  DISC 
300  DREAD(10,(Q-1)*51*1,NC13,51,E) 

310  REM  DEFINE  COMPLEX  CONSTANT  (1,01 
320  CPAK(1,0,D2) 

430  REM  ADAPTER  CORRECT  OUT  DATA 

440  FOR  1-1  TO  N 

460  CSUB (LC I  I  ,SC  I  I  ,D1 1 

480  CPAK (REA (LC I  1 1  , - I  MG ( L C I  3 1  , N2 1 

500  CMPYISCI3 ,N2,N1) 

520  CSUB (D2 ,N1 ,N1 ) 

540  CDI V (Ml ,D1 ,N3) 

553  CPAK (HAG (N3) ,0 ,N4) 

556  CDI V <N3 ,N4 ,N3) 

560  CSUB (MC I  3  ,LC I  3  ,N1 1 
580  CMP Y (MC I  3  ,N2  ,D1  1 
680  CSUB (D2 ,D1  ,D1 ) 

620  CDIV (N1 ,D1 ,N1 ) 

640  CMPY (Ml  ,N3  ,MC I  3 ) 

660  NEXT  I 

9000  REM  AUTOMATICALLY  LOAD  NEXT  PROGRAM 
9990  CHAIN (" 1 1ANU82 .S" ) 


Fig.  34.  Listing  for  Program  5B 
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Program  6  (previously  discussed)  is  then  chained 
into  memory,  unchanged,  and  the  results  are  output  in 
graphic  and/or  tabular  form. 

Thus,  lines  can  now  be  compared  directly  against 
one  line  that  serves  as  a  standard,  without  need  for  a 
perfect  termination  or  sliding  load  in  the  transmission 
line  adopted  as  a  standard:  a  procedure  not  heretofore 
available. 


RESULTS 


Four  different  SMA  lines  were  used  to  evaluate  the 
techniques  presented  in  this  document. 

1.  A  Narda  Serial  #W11  3015292G1,  which  was  fitted 
with  Narda  4401  Female  connectors  on  both  ends, 
had  an  electrical  length  of  89.4  cm,  obtained 
from  Tufts'  stock.  This  cable  will  be  titled 
"our"  line. 

2.  and  3.  Two  lines  were  constructed  at  Lincoln 

Laboratory  from  Uniform  Tubes  SMA  coaxial 
cable  fitted  with  an  OSM  210-1  male  connector 
on  one  end  and  an  OSM  207-97 76SF  female  con¬ 
nector  on  the  other.  The  "short"  line  had  an 
electrical  length  of  113.9  cm  while  the  "long" 
one  was  228.1  cm  in  length. 

4.  The  fourth  cable,  also  constructed  at  Lincoln 
Laboratory,  was  unusual  in  the  sense  that  it 
was  constructed  from  a  173.6  cm  electrical 
-length  of  Precision  Tubes  141  series  air- 
articulated  3.5  mm  line  fitted  with  Solitron/ 
Microwave  2902-6057  male  connectors  on  both 
ends.  This  line  was  not  dielectric  filled  and 
was  called  the  "spline"  line. 
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a. 

A 


These  four  lines  offered  the  chance  of  testing  three 
conventional  SMA  lines,  two  of  which  were  identical  in 
construction  and  one  which  was  a  significant  departure  in 
design. 


The  first  set  of  results  to  be  displayed  will  be 
the  adapter-corrected  reflections  of  a  zero-plane  short, 
an  HP  offset  4  short,  and  an  open  circuit  referenced  to 
the  running  averaged  load  characterization  of  each  line 
measured  from  2-17  GHz.  A  female  SMA  adapter  was  used  at 
the  measurement  port  for  all  the  measurements  and  a 
standard  SMA  fixed  load  was  used  to  terminate  each  line  for 
load  characterization. 

As  can  be  seen  in  Figs.  35-38,  the  measurement  of 
the  zero-plane  short  circuit  was  virtually  identical  for 
all  the  lines,  with  a  0.03  maximum  deviation  from  unity 
reflection  coefficient  magnitude.  Some  of  this  deviation 
may  be  attributed  to  the  slight  power  dissipation  in  the 
adapter,  which  is  assumed  tc  be  dissipationless.  The 
variations  are  probably  caused  by  the  frequency  varying 
directivity  of  the  directional  coupler.  It  is  speculated 
that  a  coupler  which  is  smoother  in  directivity  over  the 
frequency  band  (or  a  broadband  SWR  bridge,  such  as  a 
Wiltron  58A50)  would  be  more  satisfactory  in  this  respect. 
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Figure  39  shows  a  plot  of  the  coupler  directivity 
for  the  HP  8743  test  set.  This  was  measured  by  placing 
an  APC-7  sliding  load  at  the  reflection  measuring  port, 
and  by  placing  the  ANA  in  manual  mode.  A  time-exposed 
photograph  of  the  HP  180A  rectangular  display  was  made  as 
the  sliding  load  was  moved  from  stop  to  stop  several  times 
for  each  frequency  band.  The  bisection  of  the  envelope 
was  taken  to  be  a  measure  of  the  directivity.  As  shown, 
the  directivity  does  not  everywhere  safely  meet  the  as¬ 
sumption  of  being  varying  slowly  with  frequency.  This 
deficiency  could  produce  a  small  error  in  magnitude 
measurements  even  for  a  short  circuit. 

At  this  point,  the  system  stability  has  been  shown 
to  be  sufficient  for  adapter  correction.  Now,  the  aver¬ 
aged  load  and  line  and  zero-plane  short  circuit  will  be 
used  as  standards  while  two  conventional  calibration 
standards  (offset  4  short  circuit  and  open  circuit)  will 
be  used  as  DUT's  in  order  to  show  the  suitability  of  the 
complete  calibration  procedure. 

The  magnitude  and  phase  plots  of  the  offset  4  shorts 
for  each  of  the  lines  is  presented  in  Figs.  40-47.  They 
are  unremarkable  in  that  they  are  essentially  the  same  for 
each  of  the  lines. 
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Fig.  41.  Reflection  of  offset  4  short 

referenced  to  "our"  line  (Ang) . 
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Reflection  of  offset  4  short 
referenced  to  "short"  line  (Mag) . 
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Fig.  45.  Reflection  of  offset  4  short 

referenced  to  "long"  line  (Ang) . 
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It  is  important  to  note  that  the  lack  of  difference 
between  the  lines  for  the  zero-plane  and  offset  4  short 
circuits  is  due  to  the  fact  that  these  are  high  reflection 
standards  with  a  well  defined  plane  of  measurement.  Since 
zero-reflection  correction  provided  by  the  calibration  is 
so  small,  it  becomes  almost  negligible  in  measurements  of 
high  reflection  standards.  Therefore,  the  characteristics 
of  the  different  lines  used  for  the  running  averaged  load 
characterization  have  little  affect  on  the  measurements 
of  these  standards.  However,  as  will  be  shown  next,  the 
line  characteristic  will  have  a  dramatic  effect  on  devices, 
such  as  open  circuits,  where,  even  though  they  are  high- 
reflection  devices,  they  exhibit  a  frequency-varying  phase 
which  is  dependent  upon  the  nature  of  the  interconnecting 
transmission  line. 

It  will  also  be  seen  that  the  characteristics  of 
lines  can  be  significantly  different  even  among  those  of 
the  same  type  and  can  be  readily  observed  by  employing  the 
direct  line  comparison  technique  previously  described. 

Figures  48-55  show  the  corrected  measurement  of 
reflection  magnitude  and  phase  of  an  open  circuit  for  each 
of  the  lines. 
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Fig.  48.  Reflection  of  open  circuit 

referenced  to  "our"  line  (Mag) . 
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Fig.  49.  Reflection  of  open  circuit 

referenced  to  "our"  line  (Ang) 
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Fig.  53.  Reflection  of  open  circuit 

referenced  to  "long"  line  (Ang) . 
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Fig.  55.  Reflection  of  open  circuit 

referenced  to  "spline"  line  (Ang) . 


The  magnitude  plots  are  virtually  the  same,  but  as 
would  be  expected,  the  phase  plots  are  quite  different. 

Even  the  "short"  and  "long"  line  exhibited  different  phase 
characteristics  and  they  are  of  the  same  construction 
except  for  length.  There  is  also  no  apparent  correlation 
between  the  phase  variations  of  the  open  circuits  of  the 
same  sex  either.  The  phase  characteristics  of  the  SMA 
open  circuit  are,  therefore,  more  dependent  upon  the 
interconnecting  line  characteristics  than  on  the  general 
construction  of  the  line  or  the  connectors.  This  evidence 
supports  the  earlier  claim  that  it  would  be  practically 
impossible  to  characterize  an  SMA  open  circuit  without 
specifying  the  transmission  line  defining  zero  reflection. 

The  final  set  of  results  contain  the  reflection 
plots  for  the  direct  comparison  of  each  transmission  line 
against  each  of  the  others  for  14-16  GHz  region. 

Figures  56-57  show  the  magnitude  and  phase  plots  of  the 
comparison  "our"  line  versus  "our"  line.  Indeed,  the 
magnitude  and  phase  are  zero  as  would  be  expected  when 
comparing  two  identical  lines.  Figures  58-81  show  the  real 
and  imaginary  components  of  the  reflection  for  each  of  the 
other  comparisons. 

Three  noteworthy  observations  can  be  made.  First, 
two  lines  of  the  same  construction,  but  different  lengths, 
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Comparison  of  "our"  line  to 
"our"  line  (Mag) . 
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Fig.  60.  Comparison  of  "our"  line  to 
"short"  )ine  (Real). 
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Fig.  66.  Comparison  of  "spline"  line  to 
"our"  line  (Real) . 
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Fig.  67.  Comparison  of  "spline"  line  to 
"our"  line  (Imag) . 
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6R.  Comparison  of  "our"  line  to 
"spline"  line  (Real) . 
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Comparison  of  "long”  line  to 
"short"  line  (Real) . 
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Fig.  71.  Comparison  of  "long"  line  to 
"short"  line  (Imag) . 
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Fig.  72.  Comparison  of  "short"  line  to 
"long"  line  (Real) . 
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Fig.  73.  Comparison  of  "short"  line  to 
"long"  line  (Imag). 
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74.  Comparison  of  "spline"  line  to 
"short"  line  (Real) . 
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76.  Comparison  of  "short"  line  to 
"spline"  line  (Real) . 
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Fig.  77.  Comparison  of  "short"  line  to 
"spline"  line  (Imag). 
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would  be  expected  to  match  well.  Looking  at  the  results 
of  the  "long"  vs.  "short"  line,  one  can  see  that  they  do 
match  quite  well,  but  they  are  also  not  identical. 

Second,  when  comparing  the  "long"  versus  "short" 
and  then  the  "short"  versus  "long"  (or  any  other  analagous 
sets  of  data) .  One  can  see  that  real  and  imaginary  parts 
for  one  comparison  are  reflected  from  zero  in  the  opposite 
comparison,  as  would  be  expected. 

Third,  it  is  clear  that  the  characteristics  of  the 
"spline"  line  are  very  much  different  from  those  of  the 
other  three  lines.  On  the  contrary,  the  characteristics 
of  the  other  three  lines  do  not  differ  widely.  Apparently 
the  "spline"  line  does  not  have  a  50  +  jO  ohm  character¬ 
istic  impedance. 
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v  , 


EXTENSION  OF  ADAPTER  CORRECTION  TO  TWO-PORT  DEVICE  DE¬ 
EMBEDDING 

Up  to  this  point,  discussion  has  been  limited  to 
1-port  reflection  measurements.  A  large  proportion  of 
microwave  devices  are  two-port  in  construction  and  cannot 
be  fully  characterized  by  reflection  measurements  only. 
Therefore,  a  method  of  extending  the  adapter  correction  to 
two-port  measurements  needed  to  be  developed.  One  method 
will  now  be  described. 

In  connection  systems  such  as  SMA  there  are  two 
types  of  two-port  devices:  insertable,  which  have  one 
female  and  one  male  connector,  and  noninsertable ,  which 
have  the  same  sex  connectors  at  both  ports,  be  they  male 
or  female.  As  a  consequence,  any  adapter  kit  used  for 
general  two-port  measurements  must  include  two  adapters 
for  each  sex  and  the  same  sex  adapters  cannot  be  assumed 
to  be  identical.  Therefore,  a  method  of  measuring  the 
scattering  matrices  of  these  four  adapters  in  order  to 
de-embed  the  scattering  matrix  for  the  DUT  will  now  be 
presented. 

Let  A  and  B  represent  the  scattering  matrices  for 
the  adapters  of  one  sex  and  C  and  D  of  the  other.  It  will 
be  reasonably  assumed  that  all  adapters  are  passive  and 
reciprocal,  that  a  proper  APC-7  calibration  has  been  ac- 
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complished  and  that  all  measurements  are  corrected  for 
system  errors  against  this  calibration.  (Note  that  the 
DUT  need  not  be  reciprocal.)  Also,  discussion  will  assume 
that  SMA  adapters  are  being  used,  but  the  method  is  di¬ 
rectly  applicable  to  any  type  of  adapter. 

Letting  the  APC-7  end  of  one  sex  adapter,  say 
female,  be  port  1,  the  SMA  end  be  port  2,  the  SMA  end  of 
the  male  adapter  be  port  3,  and,  finally,  the  APC-7  end  of 
the  male  adapter  be  port  4,  one  obtains  the  following 
f lowgraphs : 


Female  Male 


Fig.  82.  Flowgraphs  for  one  set  of  four 
adapters . 


157 


All'  Bll'  C44'  and  D44  can  be  found  directly  by 
measuring  the  reflection  of  each  adapter  when  terminated 
by  a  matched  load  (fixed,  sliding,  or,  as  previously 
discussed,  the  Running  Averaged  load  and  line)  on  the  SMA 
end. 


Reflection  measurements  are  next  made  on  each 
adapter  with  a  zero-plane  short  circuit  attached  to  the 
SMA  end.  The  following  equations  result  where  F^,  T2>  F^, 
and  P4  represent  the  measurements  on  the  A,  B,  C,  and  D 
adapters,  respectively. 


r.  =  a.  .  + 


=  B,,  + 


r,  =  c.,  + 


r .  =  d  . .  + 


(A2l) 

2  (-1) 

1  + 

A22 

(B21> 

2  (-1) 

i  + 

<M 

CM 

CO 

<C43> 

2  (-D 

1  + 

C33 

(D43) 

2  (-1) 

1  + 

D33 

(32) 


(33) 


(34) 


(35) 
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2  2  2 

By  rearranging  Eqs.  32-35  and  solving  for  A^,  B21'  ^43 
2 

and  one  obtains 


(An  -  Ti)  (1  +  A22) 


(36) 


B 


2 

21 


r2m  ♦  b22) 


(37) 


C43  *  <C44  -  r3l(1  +  C33>  (38) 

D43  *  <D44  '  r4’  (1  +  “SS1  1391 


Now  attaching  the  SMA  ends  of  A  and  C  and  of  B  and  D 
together,  one  makes  reflection  measurements  on  each  end  of 
each  adapter  pair  with  the  other  end  terminated  with  a 
matched  load,  obtaining 


r 


5 


21  33 


1  -  A 


22  33 


(40) 


■3* 
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c2  A 

r  =  43  22 

l6  1  -  C. 


'22  33 


(41) 


r,  = 


B21  D33 
1  "  B22  D33 


(42) 


r„  » 


P43  B22 
1  "  B22  D33 


(43) 


where  r,.,  rg ,  r7  and  Tg  are  the  measured  reflection.  Re¬ 
arranging  Eqs.  40-43  yields 


21 


21 


'43 


'43 


r5(l 

"  A22 

C33' 

C33 

r7  ( 1 

B22 

D33> 

D33 

V1 

A22 

C33} 

A2  2 

r8(1 

B22 

°33) 

(44) 


(45) 


(46) 


22 


(47) 
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Appropriately  equating  Eqs.  36-39  to  Eqs.  44-47 
rearranging  terms  yields 


r5 

.  C33(1  + 

a22) 

=  A11 

“  F1 

1  "  A22 

C33 

= 

r7 

_  °33(1  + 

B22  ^ 

-  r2 

1  "  B22 

°33 

r6 

A22^  + 

C33) 

C44 

-  F  3 

1  "  A22 

C33 

r8 

_  B22<1  + 

°33) 

Q 

1 

-  r4 

1  "  B22 

°33 

where  X^,  X2 ,  x^,  and  X^  are  calculable  as  shown 
for  2  in  Eq.  48,  D33  in  Eq.  49,  A22  in  Eq.  50, 
in  Eq.  51,  one  finds 

xi 

C33  '  X  ♦  A22  *  Xj  A22 


and 


(48) 

(49) 

(50) 

(51) 

Solving 
and  B22 


(52) 
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X, 


'33 


!  +  B22  +  X2  B22 


(53) 


22 


1  +  C33  +  X3  S3 


(54) 


22 


1  +  D 


33 


+  X, 


33 


(55) 


Substituting  in  Eq.  52  into  Eq.  54  and  ^  in  Eq.  53 

into  Eq.  55,  we  find 


A 


22 


(56) 


22 


4 


(57) 


Substitution  of  A22  and  B22  intc  Eqs.  52  and  54  permits 
calculation  of  and 

Since  two  possible  square  roots  exist  for  any 
Comdex  number,  each  separated  by  180°  from  the  other  in 


-M 

I 


the  complex  plane,  it  first  appeared  that  a  method  of 

choosing  the  correct  root  would  be  necessary.  However,  it 

was  observed  that  since  the  adapters  are  always  used  in 

pairs,  and  since  the  roots  differ  only  in  sign,  it  is  not 

necessary  to  solve  for  the  correct  root.  All  that  needs 

2  2''  2 

to  be  known  is  that  the  roots  of  A^.,  ,  C*!  - ,  and  _ 

21  21  43  43 

are  all  correct  or  that  all  are  incorrect  because  the 

products  A2lC43'  A21D43'  etc*  rePr©senting  the  trans¬ 
mission  paths  through  the  adapter  pairs  will  be  equal  in 
either  case.  This  observation  greatly  simplifies  the  task 
of  finding  the  roots. 

One  method  is  to  measure  the  transmission  through 

adapter  pairs  combined  in  three  different  ways.  For 

example,  let  Tj.  =  A21C43,  T2  =  B21D43,  and  T3  =  B2]C43. 

Then  solve  for  one  root  of  A_,  ,  B_,,  C.,,  and  D .  (A  ,  B  , 

21'  21  43'  43  r'  r ' 

Cr,  and  D^,  respectively). 

If  angle  (Arcr)  ^  angle  (T^) ,  then  reverse  the  sign 
of  either  or  C^,.  Otherwise,  accept  the  signs  as  they 
are.  Secondly,  if  angle  (B^C^.)  ^  angle  (T3)  ,  then  reverse 
the  sign  of  B^;  otherwise  accept  it.  Thirdly,  if  angle 
(B^D^.)  ^  angle  (T2)  ,  reverse  the  sign  of  D^;  otherwise, 


accept  that  one. 
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At  this  point,  either  the  signs  of  Ar,  Br,  C^,  or 
are  all  correct  or  they  are  all  incorrect,  exactly  the 
information  needed  to  proceed  with  the  de-embedding  of  the 
scattering  matrix  of  a  DUT  connected  between  any  pair  of 
adapters.  It  is  important  to  note  that  the  numbering  of 
the  ports  has  been  chosen  for  convenience  here,  and  that 
they  can  be  represented  in  any  suitable  manner.  Note  that 
the  adapters  are  reciprocal  and  that  the  forward  and 
reverse  transmission  parameters  of  each  adapter  have  been 
freely  interchanged. 

Now  that  the  scattering  matrices  have  been  found 
for  each  adapter,  de-embedding  of  ^he  scattering  parameters 
for  the  DUT  alone  can  be  performed.  By  de-embedding  one 
adapter  at  a  time,  one  avoids  solving  a  pair  of  complex 
quadratic  equations  and  an  explicit  solution  can  be 
obtained  with  gratifying  simplicity.  This  procedure  is 
much  easier  to  understand  and  implement  than  a  previous 
technique  presented  by  Saleh  (Ref.  14) . 

For  the  sake  of  discussion,  let's  assume  that  a 
scattering  matrix  E  is  measured  for  a  DUT  in  combination 
with  adapters  A  and  C,  as  above,  even  though  any  pair  of 
adapters  can  be  used.  Treating  the  DUT  and  adapter  C  as  a 
lumped  device  with  scattering  matrix  £ ,  E  can  be  repre¬ 
sented  by  a  model  as  illustrated  below. 


Fig.  83.  Model  for  separating  adapter  A 
from  lumped  device  £. 

From  the  flowgraph  for  the  model  above  the  following 
equations  result. 


(59) 


V 
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Lastly,  it  is  necessary  to  de-embed  the  scattering  matrix 
S  for  the  DUT  from  /.  The  following  model  is  then  employed. 


Fig.  84.  Model  for  separating  adapter  C 
from  DUT. 

From  this  model  one  finds 


„  =  c  .  S21  S12  C33 
'll  -  S11  1  -  S22  c33 


(66) 


1 


S2 1  C4  3 
S22  C33 


(67) 


12  C43 


S22  C 


33 


(68) 


As  can  be  seen,  the  explicit  solution  for  scattering  matrix 
S  can  be  obtained  from  these  equations  yielding 


S 


+  *22  C33  ‘  C33  C44 


(70) 


=  *21(1  ~  S22  C33} 
21  “  C43 


(71) 


S 


12 


*12  ( 1  -  S22  C33> 


'43 


(72) 


S 


11 


(73) 


Thus,  it  is  clear  that  by  using  this  method  of 
adapter  correction  and  de-embedding  for  two-port  devices, 
the  scattering  matrix  for  any  DUT  can  be  obtained  without 
iteratively  solving  simultaneous  complex  quadratic  e- 
quations  and  only  with  reference  to  zero-plane  short- 
circuit  and  matched  load  standards. 


CONCLUSION 
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What  have  been  presented  here  are  new  and  inno¬ 
vative  techniques  for  calibrating  Automatic  Microwave 
Network  Analyzers  (e.g.  HP  8542B)  to  make  complex  re¬ 
flection  measurements.  These  techniques  offer  significant 
advantages  particularly  when  making  measurements  using 
miniature  connector  systems. 

It  has  been  shown  that,  by  using  the  experimentally 
derived  characterization  of  an  APC-7  open  circuit  as  a 
high  reflection  standard,  band  edge  discontinuities  can  be 
eliminated  with  no  loss  of  measurement  accuracy.  Use  of 
multiple  offset  shorts  is  no  longer  necessary  resulting  in 
fewer  connection  requirements.  This  decrease  translates 
into  less  wear  of  measurement  port  connectors  or,  more 
basically,  lower  maintenance  costs. 

The  adapter  correction  based  on  low  dissipative 
losses  permits  measurements  in  other  connection  and  trans¬ 
mission  formats  without  the  need  for  a  full  calibration 
kit  or  process  in  the  adapted  format.  This  greatly  reduces 
operating  costs,  requiring  fewer  standards  and  connections 
especially  in  the  delicate  and  easily  damaged  SMA  trans¬ 
mission  system. 

Finally,  the  use  of  a  known  length  of  terminated 
transmission  line  and  the  Double  Running  Average,  instead 
of  a  sliding  load  or  fixed  termination,  to  characterize  a 
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perfectly  matched  load  offers  a  means  of  establishing  a 
low-reflection  standard  truly  relevant  to  the  measure¬ 
ment  medium  when  using  dielectric  filled  lines  for  inter¬ 
connection.  This  technique  eliminates  the  need  for  air 
sliding  loads  in  the  adapted  formats  (particularly  the 
delicate  SMA) ,  excludes  the  possibility  of  encountering 
defective  calibration  due  to  divergences  of  the  circle 
fitting  algorithm  and  enables  direct  comparisons  of  differ- 
erent  transmission  lines. 

Thus,  at  the  relatively  insignificant  cost  of 
increased  computer  memory  and  slightly  longer  calibration 
times,  important  cutbacks  in  maintenance  and  equipment 
costs  and  greatly  decreased  operator  assistance  during 
measurements  can  be  realized  for  existing  types  of 
measurements.  Measurements  of  exotic  transmission  lines, 
such  as  microstrip,  slotline,  and  dielectric  image  guide, 
and  of  any  transmission  formats  developed  in  the  future, 
are  not  only  possible,  but  are  easily  accomplished  without 
the  need  for  perfect  or  sliding  loads. 
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APPENDIX 

This  appendix  contains  the  listings  of  the  special 
software  mentioned  in  the  chapter  entitled  Frequency 
Characterization  of  APC-7  Open-Circuit  Reflection.  These 
programs  used  the  HP  error-correction  algorithms,  but  they 
allowed  specification  of  offset  lengths  for  the  zero-and 
X/4-planes. 

Also  contained  in  this  appendix  are  sectioned  views 
of  two  typical  SMA  connectors. 


i 


pace  l 

1  PEM  THIS  PROGRAM  DEVELOPED  OH  2/20/80 .  IT  USES  COMPLEX 

2  PEM  AVERAGED  SLID  SLIDING  LOAD  AND  OFPSET  REFERENCE  PLANE 

3  PEM  SHORTS  AS  REFERENCES  TO*MEASURE  OPEN  CIRCUIT  REFLECTION 

4  REM  PEV  3x6/31  REM  3/6/81  1  i  A N U 4 1  .  S 

20  COM  K C251 3  .BC251 1  .RC231 3  , AC231 3  .UC231 3  .TC231  3  .FC331  3  ,NC 103 
40  DSPLAY  "FREO(MHZ)  -  S T A R T  ,  S T OP  ,  S T E P  “  1 
50  BELL 

60  INPUT  F1.F2.S 

70  REM  CALCULATE  NUMBER  OF  FREQUENCIES  OF  INTEREST 
80  LET  N« 1 ♦ < F2-F 1 ) x$ 

33  REM  INITIALIZE  RF  MEASUREMENT  EQUIPMENT 
30  FCALF(Fl) 


100 

BCNTl  (FI  ) 

I 

1  20 

33EL1 (11) 

1  40 

WAIT  (30) 

150 

REM  CONNECT  APC-7  REFERENCE  PLANE  SHORT  AND  MAKE 

151 

PEM  MEASUREMENTS 

- 

1  60 

DSPLAY  "CONNECT  REFERENCE  PLANE  SHORT” 

1  30 

PAUSE 

200 

FOP  1-1  TO  N 

210 

REM  CALCULATE  AND  STORE  FREQUENCIES  OF  INTEREST 

220 

let  rm.Fi+ii-ius 

225 

PEM  INIALIZE  TO  ZERO  REAL  AND  IMAGINARY  PARTS  OF 

SLIDING 

226 

REM  LOAD  COMPLEX  SUM 

2  30 

LET  UCI3-0 

235 

LET  Tcn-o 

237 

REM  MAKE  MEASUREMENTS  ON  REFERENCE  PLANE  SHORT 

| 

240 

FPEQ2 iFtlll 

1 

260 

MEAS 1 ( 1 50 , X , Y ) 

1 

230 

CP AK (X , Y , AC  I 3 ) 

1 

300 

NEXT  I 

i 

1 

3  1  0 

REM  CONNECT  SLIDING  LOAD  AND  MAKE  6  MEASUREMENTS 

PT  EPCH  1 

31  1 

PEM  FREQUENCY  OF  INTEREST  AND  ACCUMULATE  SUMS  OF 

RE  PL  PND 

312 

REM  IMAGINARY  PARTS 

320 

DSPLAY  "CONNECT  SLIDING  LOAD” 

330 

PAUSE 

340 

FOR  J-l  TO  6 

350 

IF  J-l  GOTO  400 

. 

355 

REM  SLIDE  SLIDING  LOAD  LAMBDA/ 1 2  FOR  PRIMARY  FREQUENCY 

360 

DSPLAY  " SL I DE" 

1 

380 

PAUSE 

i 

4  00 

FOP  I -1  TO  N 

1 

420 

FPE02 ( F  C I  I  ) 

i 

4  4  0 

MEA 31(130, X,Y) 

| 

4  6  0 

let  ucn-uciJ+x 

I 

i 

4  30 

LET  t  C I  3  - T  C I  3 ♦ Y 

5O0 

NEXT  I 

520 

NEXT  J 

5  30 

PEN  PE  P  F  OPM  COMPLEX  AVERAGE 

540 

FOR  I • 1  TO  N 

i 

5  6  0 

C  P  AK <  UC  I  3  ,T  C  I  3  ,  BC  I  3  ) 

530 

C  P  A  K  (  6 . 0  ,  X  ) 

6  00 

C  D  I  V  (  B  C  I  3  ,  X  .  B  C  I  I  1 

6  2  0 

LET  UC  I  3 -PEA  1 B  C  I  3  ) 

LET  T  C  I  3  - IM  G i B  C  I  3  i 

ri  1*1 

NEXT  I 

2  5  0  0 

-'Ell  :  T  OF  E  LOCAL  VARIABLES  IN  COMMON 

£  5  *  0 

LET  1 1 ; 1  3 . F  1 

i.  *  3  0 

LET  II C  2  3  -  F  2 

i  4  I  0 

LET  NC63-S 

I  4  7 

LET  n  C  7 3 -N 

L.  4  0 

/Ell  AUTOMATICALLY  LOAD  ''EXT  PPOGPAM 

~j  t 

h  A  I  N  1  '  1  :  AINJ4  2  .  :  "  > 

Fig.  Proqram  1  for  open  circuit  phase 

measurement. 
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PAGE  1 

1  REM  THIS  PROGRAM  DEVELOPED  ON  6/3/79.  IT  IS  A  CONTINUATION 

2  REM  OF  1  : A  N  U  4 1  . C  FOR  MEASURING  THE  REFLECTION  OF  AN  OPEN 

3  REM  CIRCUIT 

4  REM  REV  3/6/81  REM  3'6/81  liANU42.S 

80  COM  Ct251.33.Lt2513.OC2513,SC2513,UC2313,TC2513.FC251),NC103 

30  REM  LOAD  LOCAL  VARIABLES  FROM  COMMON 

40  LET  I  =  N  C  6  3 

80  LET  F l -NC  1  3 

1O0  LET  F2-NC  2 3 

150  LET  N.Niri 

180  REM  INITIALIZE  RF  MEASUREMENT  EQUIPMENT 
170  FCALFIFP 
1  SO  BCNT  1  i F 1  ) 

200  SSEL  1  ill) 

220  WAIT  (50) 

230  REM  MAKE  MEASUREMENTS  ON  OFFSET  SHORT 
240  DSPLAV  "CONNECT  OFFSET  SHORT" 

280  BELL 
230  PAUSE 
320  FOR  I-l  TO  N 
340  FREQ2 ( F  t I  3 ) 

3.60  MEAS1  <  150  ,)<  ,V) 

330  CPAK (X ,Y ,0C I  3  I 
400  NEXT  I 

410  REM  READ  IN  DIFFERENCE  IN  LENGTHS  BETWEEN  REFERENCE 

411  REM  PLANE  SHORT  AND  OFFSET  SHORT 

420  DSPLAV  "WHAT  IS  OFFSET  SHORT  POSITION  (CM)  COMPARED  TO" 

42!  DSPLAV  "REFERENCE  PLANE  SHORT") 

440  BELL 
460  INPUT  C  9 

470  REM  CONVERT  LENGTH  TO  RADIANS  PER  MHZ 
430  LET  C9»-C9«2*2 .0965E-O4 

490  REM  INITIALIZE  (1,0)  AND  (-1,0)  CONSTANTS 
500  CPAK(1,0,D2) 

520  CPAK  (-1  ,0,1)31 
530  FOP  1.1  TO  N 

5  90  REM  CORRECT  OFFSET  SHORT  PHASOR  TO  LAMBDA/1 4  POSITION 
591  REM  COMPARED  TO  REFERENCE  PLANE  SHORT 

6O0  p:FT(0,FCI3,C9,D3,P1) 

660  REM  CALCULATE  ERROR  CORRECTION  COEFFICIENTS 

661  REM  C ( I  ,  1  ) »  E  0  0  ,  C ( I  , 2 ) - E  0 1  ,  C(I, 31-Ell 

6  8  0  C P  A  K  i U  C  I  3  , T  C I  3  , C  C  I  ,  l  I  ) 

7  0  0  u : U  E  <  S  t  I  3  , C I  I  ,  1 3  ,  N l ) 

720  c:UB'OCI3,CCI,l3,Dl) 

740  CDIViMt.Dl.Nl) 

7  6  0  C  M  P  V l P 1  , N 1  ,  Q  5 ) 

3  0  0  C  A  D  D i 0  5 , D 2  .  N 1  ) 

3  2  0  C S  U  B i P 1  , Q 5  ,  D l  ) 

3-0  C  D I  V i N 1  ,  D 1  , C t I  ,  3  3 ) 

3  6  0  C S  U l ' C  C  I  .  1  3  .  2  C I  3  ,  N 1 ) 

330  CADDiD2,CCI.3),Dl) 

9  00  C A D  D i  D  2  .  C  C  I  ,  3  3  ,  D  1  I 
9 2  0  C H  F  v  I N 1  . D 1  . C  C  I  .  2  3  I 
9  4  0  :  u  e  i  •_  c  1 3 . :  c  1 3 .  d  i  ■ 

■9  6  0  NE.-.T  1 

I860  C  H  A  I  M  '  “  1  i  HINJ4  3  .  I  "  1 


Fig.  86.  Program  2  for  open  circuit  phase 
measurement. 
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1  'EM  THIS  P POGROM  DEVELOPED  ON  1 1  ' 1 5  ✓  7  8  .  IT  IS  A  CONTINUATION 

2  REM  OF  1:ANU42.S  FOR  MEASURING  THE  REFLECTION  OF  AN  iiREH 

3  REM  CIRCUIT 

4  REM  REV  3 '6  '  81  REM  3 ' 6 ' 8  1  1  :  fl  N  U  4  3  .  $ 

20  COM  CC231,33,KC2313,RC2313,SC2313.UC2313,TC2313PFC2513,NC103 
SO  REM  LOAD  LOCAL  VARIABLES  FROM  COMMON 
40  LET  3.NC63 
43  LET  Fl.Nt  I  ] 

33  LET  F 2 • N  C  2  3 
7 0  LET  N »  N  t  7  I 

SO  REM  INITIALIZE  RF  MEASUREMENT  EQUIPMENT 
?  0  F  C  A  L  F  t  F  1  ) 

100  eCNTl (Ft ) 

12  0  3  S  E  L  1  <  1  1  i 
14  0  Id  A  I  T  (30) 

130  REM  CONNECT  OPEN  AND  MAKE  MEASUREMENTS 
160  DSPLAY  "CONNECT  OPEN" 

180  PAUSE 

200  FOR  I»1  TO  N 

220  FREQ2 (FC  I  1  ) 

240  ME AS  1  ( 130 ,X,Y) 

260  CPAK ( X , Y  ,  SC  I J ) 

280  NEXT  I 

280  PEM  READ  IN  LENGTH  OF  SHORT  USED  FOR  REFERENCE  PLANE 
300  DSPLAY  "WHAT  IS  POSITION  OF  REFERENCE  PLANE  (CM)  **  : 

320  BELL 
340  INPUT  Ml 

330  PEM  CONVERT  LENGTH  TO  RADIANS  PER  MH2 
360  LET  M1»-M1»242 .0963E-04 
2100  REM  INITIALIZE  (1.0)  CONSTANT 
2200  CPAK (1 .0 .D2) 

2210  REM  CORRECT  OPEN  CIRCUIT  MEASUREMENTS  UITH  ERROR 

2211  REM  CORRECTION  COEFFICIENTS 
2220  FOR  1.1  TO  N 

2240  CSUB  CSC  I  3  ,CE  I  ,  1  J  ,D1  ) 

2260  CDIVCCCI  ,23  ,DI  ,01) 

2280  C  A  D  D ( D 1  ,Ct I  ,3]  ,D1 > 

2300  CDIVCD2.D1  ,SCI3) 

2310  PEM  ROTATE  CORRECTED  OPEN  CIRCUIT  MEASUREMENT  TO  O-PLANE 

2311  PEM  REFERENCE 

2320  PSFT (0  .FC  I  3  ,M1  ,3C I  3  ,SC  I  3  I 
2340  NEXT  I 

2  3  00  CHAINi  Mi  IANU44.S”) 


Fig.  87.  Program  3  for  open  circuit  phase 
measurement. 
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PAGE  1 

1  REM  THIS  PROGRAM  DEVELOPED  OH  10'11'79  PHD  IT 

2  REM  OUTPUTS  THE  DPTP  GENERATED  BY  ANU41.S,  ANU42.S,  AMD 

3  REM  ANU43.S  TO  THE  VIDEO  PLOTTER 

4  REM  REV  5/-l'81  REM  3-'20-'31  J:PNU44.S 

20  COM  AC251 I  ,BC231  I  ,01231  ,31  ,KC251 3  ,Rt231  3  ,MC231 J 

30  COM  UC251  3  ,TC231  3  ,F[23I  3  ,NC  103 

35  REM  RECALL  LOCAL  VARIABLES  PROM  COMMON 

48  LET  S  »  N  C  6  3 

120  LET  F5=NC13 

140  LET  F6-N  C  2  3 

150  LET  N.NC73 

130  REM  INITIALIZE  GRAPHICS  DISPLAY  AND  BUFFER 
2  0  0  B  U  F i  1 0  0  i 
220  CLEAR  1 0  > 

230  DSPLPY  “UHAT  IS  B1  AND  Tl") 

235  INPUT  Bl.Tl 

260  LET  3 1  - i  T 1 -B  1  ) y  10 

270  REM  CLEAR  CRT,  SCALE,  AND  PLOT  REAL  PART  OF  REFLECTION 

230  clear  1 0 i 

660  SCALE iF3,F6,Bl,Ttl 
680  LET  F3- ( F6-F3 ) ' 1 0 
700  SAXES (  F  3 ,31  ) 

710  LABEL  VTAB ( 13) ,TAB (01 , "REFL” 

720  LABEL  VTPB(16) ,TAB(0) ."REAL" 

740  LAEEL  VTPB ( 3  1  )  , TAB (28)  , " FREQ ( MHZ)  ■ 

760  BL0CK(B2) 

780  FOP  I»1  TO  N 

800  PLOT ( F [ I  3  ,REA (MC I  3 )  ,2 ) 

820  NEXT  I 
920  BELL 
940  PAUSE 

930  REM  CLEAR  REAL  PLOT  AND  PLOT  IMAGINARY  PART 
960  CLEAR ( 0 ) 

970  SAXES (F3  ,S1 ) 

973  LABEL  V T A B ( 1 3 )  ,  T A B ( 0 )  ,  " RE FL " 

930  LABEL  V  T  A  B  (  1  6  3  ,  T  A  B  ( 0 )  ,"IMPG" 

1000  LABEL  VTAB ( 31 ) ,TAB (23) , “FREQ CMH2) " 

1020  BLOCK ( B 2  J 

10  4  0  FOR  I  a  1  TO  N 

I860  PLOT (Ft  I  3  , IMG (MC I  3 )  ,2) 

1080  NEXT  I 
1100  PAUSE 
1128  C L  E  A P ( 0 ) 

1140  D  SPLAY  "1.1  HAT  IS  B2  AND  T2"I 
1160  INPUT  B2.T2 

1170  PEM  CLEAR  DISPLAY  AND  PLOT  REFLECTION  MAGNITUDE 


113  0 

CLEAR ( 0  j 

1200 

LET  S2«(T2-B2)'10 

12  10 

SC«LE<F3#F*fB2 

,T21 

122  0 

:AXC?(F3,32) 

1240 

LhBEL  VTAB i 15) 

,  T  A  B  (  0  )  ,  "  R  E  F  L  " 

12  60 

LABEL  VTAB  f 161 

,  T  A  B  (  0  )  MAC*' 

1  2  3  0 

L  A  t  £  L  V  T  A  B  <  3  1  ■ 

,  T  A  B  (  2  3  )  ,  "  FPEO(MHZ)" 

1  3  0  0 

B L  0 C V  'E  2 ' 

13  2  0 

F  0  F  1*1  T  0  N 

1  3  40 

c'  L  0  T  »  F  C  I  3  .MAC  i  M  C  I  ]  >  ,2' 

1  3  6  0 

r*  e : :  T  I 

1  :  0 

F  A  IJ  '  Z 

Fig.  88.  Page  1  of  Program  4  for  open  circuit 
phase  measurement. 
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Sectioned  view  of 
OSM  207-9776SF 
Female  SMA 
connector . 


Sectioned  view  of 
Narda  4401  Female 
SMA  connector. 
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